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ABSTRACT 


An  experimental  investigation  of  the  effect  of  variation 
in  flow  angle  on  convection  heat  transfer  and  frictional 
resistance  in  a  tube  bank  i^  undertaken  ior  three  systematically 
varied  geometries.   Data  taken  by  previous  investigators  tor 
a  fourth  geometry  is  included  in  the  analysis. 

Design-oriented  curves  and  equations  derived  from  analysis 
of  the  results  are  presented.   The  curves  give  the  value  of  an 
inclination  factor  which  may  be  used  as  an  additional  factor 
in  the  modified  Grimison  equation. 

Similar  design  oriented  curves  are  presented  which  give  an 
inclination  factor  which  is  in  effect  the  ratio  of  friction 
factor  at  any  angle  to  that  for  cross  flow.   Equations  for 
the  friction  factor  at  ninety  degrees  are  derived. 

Extensive  use  is  made  of  electronic  digital  computers  for 
data  reduction  and  analysis. 


NOMENCLATURE  AND  DEFINITIONS 

Program        Engineering  Definition 

Listing  Symbol 

A  A         The  inclination  factor  equation  inter- 

cept in  the  convection  heat  transfer 
regression  analysis. 

AK  K         Flow  coefficient  for  a  sharp-edged  thin 

plate  orifice. 

ANRE  N  Reynolds  number  through  the  thin  plate 

orifice. 

ANNU  N  y        Experimental  Nusselt  Number,  Cross  flow. 

ARGFAC  F        Arrangement  factor  in  the  modified 

Grimison  equation. 

B  B        Coefficient  of  Reynolds  Number  in  the 

convection  heat  transfer  regression 
analysis. 

BK  K         Flow  coefficient  for  sharp-edged  thin 

plate  orifice  calculated  by  data 
reduction  program. 

BNRE  Nr,  „        Reynolds  Number  in  the  tube  bank  where 

RE 

transport  data  is  a  function  of  TF 

BNNU  N_T        Experimental  Nusselt  Number,  c-ny  angle 

C  -        Quotient  of  AK  and  BK  taken  so  as  to 

always  r>e  less  than  unity. 


Program 
Lifting 

CNRE 


CNNU 

CONST 

CPD 

CPH 

CPW 

DHF 

DIAFAC 

DNRE 
DP 

ENRE 


Engineering 
Symbol 


N 


N 


RE 


NU 


K 


P 


Definition 

Reynolds  Number  in  the  tube  bank  where 
transport  data  is  a  function  of  TFF. 
Nusselt  Number  calculated  by   Colburn 
equation. 

The  coefficient  of  Reynolds  Number  in 
the  frictional  resistance  analysis. 

Specific  heat  of  dry  air  at  constant 
pressure. 

Specific  heat  of  water  at  constant 
pressure. 

Specific  heat  of  wet  air  at  constant 
pressure. 

Pressure  drop  across  the  thin  plate 
orifice. 

Two  tenths  power  of  the  ratio  of 
hydraulic  diameter  to  tube  diameter. 

Floating  point  mode  equivalent  of  NRE. 

Pressure  drop  across  the  tube  bank  in 
inches  of  water. 


NR„         Reynolds  Numbers  used  in  computation 
of  friction  factor  from  frictional 
resistance  analysis. 


0.20 
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Program 
Listing 

EXPO 


FRICT 


Engineering 
Symbol 


FRIFAC 


IDO 

IRCN,  RNCD 

Equation 

KASES 


LAF 

LANGLE 

LGF 


M 


Definition 

Exponent  &f  Reynolds  number  in  the 
frictional  resistance  analysis 

Friction  factor  computed  from 
parameter  developed  by  frictional 
resistance  analysis. 
Friction  factor  computed  in  data 
reduction  program. 

Convection  heat  transfer  coefficient. 

Diameter  of  the  thin  plate  orifice 

Run  code  number  for  identification. 

Thermal  conductivity  of  air 

Maximum  number  of  data  sets  in  any 
run  of) the  data  reduction  or  frictional 
resistance  analysis  programs,  index 
of  controlling  DO  Loop  in  each. 

Back  pitch  indicator. 

Angle  indicator. 

Arrangement  indicator  index  of  DO  Loop 
representing  the  sequential  number 
of  a  set  of  data  being  processed  in 
data  reduction  program. 

Exponent  of  Reynolds  Number  in  convection 
heat  transfer  analysis. 


Program       Engineering         Definition 
Listing         Symbol 

NRE  -  Index  of  DO  loop  in  frictional 

resistance  analysis  for  computation  of 

frictional  resistance  from  parameters 

computed  in  the  program. 

NUMB  -  Inaex  analysis  program  of  DO  Loop 

in  frictional  resistance.   Maximum 
number  of  data  points  in  a  uata  set, 
in  the  frictional  resistance  analysis. 

p,  p.  Pressure  above  atmospheric  aheaa  of 

the  thin  plate  orifice 

PB  p.  Barometric  pressure 

RHO  /0  Air  density 


■•■ 


PHU  X  Specific  humiuity 


TAVG  T  vr         Average  or  bulk  temperatu; e  of  air  in 

the  tube  bonk. 

TE  T„  Exhaust  temperature;  temperature  of 

air  after  passing  through  the  tube 
bank. 

TF  T  Film  temperature,  mean  between  bulk 

temperature  and  tube  surface 
temperature . 

TFF  T  Film  temperature,  for  frictional 

resistance  calculation,  based  on 
Ts  -  C  (LMTD) 
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Program       Engineering        Definition 
Lifting         Symbol 

THUMB  -  Floating  point  moue  equivalent  of 

NUMB 

TT  T.,T  Tube  surface  temperature 

t   s  ^ 

V  -  Intermediate  quantity  calculated  in 

the  data  reduction  program. 

WAO  W  Weight  rate  of  flow  of  air  tnrough 

the  orifice. 

XCON  log  K  Natural  logarithm  of  CONST. 

3e 

Y  Y  Net  expansion  factor,  ratio  of  the 

flow  coefficient  of  a  gas  to  a 
liquid  at  the  same  Reynolds  Number. 

2  -  Intermediate  quantity  calculated  in 

the  data  reduction  program. 

-  '  Dynamic  viscosity  of  air. 


o 


Angle  of  inclination  of  tuoe  cank; 
cross  flow  is  ninety  degrees, 
parallel  flow  zero  aegrees. 
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DEFINITIONS  OF  COMPUTER  ORIENTED  TERMINOLOGY 


Fixed-point  variable  - 


Floating-point  variable  - 


DO  loop  - 


a  variable  written  without  a  decimal 

point  using  tne  digits  o,  1 9; 

may  be  positive  or  negative, 

a  variable  written  with  a  decimal  point 
using  the  digits  o,  1....9,  with  an 
optional  preceeding  plus  or  minus  sign. 
The  variable  value  may  contain  a 
multiplier  v.hich  is  ten  to  some  power, 
the  multiplier  being  represented  by 
the  symbol  E  and  the  power  of  ten. 

a  method  of  performing  the  same 
calculation  with  different  oata  by 
repeatedly  executing  the  statements 
following  tne  DO  statements  the  number 
of  times  indicated  by  an  inaex. 


BLOCK  DIAGRAM  SYMBOLS 


Processing  a  group  of  program  instruc- 
tions which  perform  a  processing 
function. 
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V 


Input  or  output  -  any  function  of  any 
input  or  output  device. 


Decision  or  Logic  -  The  decision  or  logic 
function  documents  points  in  the  program 
where  a  branch  to  alternate  points  is 
possible. 


Predefined  process  -  a  group  of  operations 
not  detailed  by  the  particular  program. 


Terminal  -  the  beginning  end  or  a  point 
of  interruption  of  a  program. 


U 


An  off-page  connector 


Magnetic  Tape 


Connector 
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RUN  CODE  NUMBERS 

A  system  of  coding  data  points  to  reflect  in  a  single 
multi-digit  number,  tube  geometry,  angle  of  inclination  and 
sequential  number  of  the  test  run  was  developed.   The  run  code 
number  is  a  four  digit  number;  the  first  number  indicating 
tube  geometry,  the  second  indicating  tne  angle  of  inclination, 
and  the  third  and  fourth  numbers  indicating  the  sequential 
number  of  the  run. 

Tube  geometry  codes: 
Code  Description 

1  Staggered  tubes,  ST  =  SL  =  1.5  D   (Smith-Kiss  and 

Bond-wallin  data) 

2  In-line  tubes  ST  =  SL  =  1.5  D 

3  Staggered  tubes  ST  =  l.S  D   SL  =  2.0  D 

4  In-line  tubes  ST  =  1.5  D    SL  =  2.0  D 

Angle  codes: 
Code  Description 

3  30  degree  angle  of  inclination 

4  45  degree  angle  of  inclination 

6  60  degree  angle  of  inclination 

7  75  degree  angle  of  inclination 
9  90  degree  angle  of  inclination 

Example:   1704 

tube  geometry  1,  ST  =  SL  =  1.5  D,,  staggered  tubes, 
seventy  five  degree  angle  of  inclination, 
fourth  run 
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INTRODUCTION 

While  adequate  empirically  derived  design  information  is 
available  for  heat  transfer  equipment  where  the  tuoe  banks  are 
parallel  to  or  normal  to  the  flow  of  the  fluid  being  heated  or 
cooled,  little  information  exists  for  the  case  where  the  tudes 
are  inclined  at  some  angle  to  the  fluid  flow.   The  case  of 
tubes  at  an  angle  to  fluid  flow  is  not  uncommon  particularly 
in  marine  boiler  convection  banks.   It  is  the  purpose  of  this 
thesis  to  investigate  the  effect  of  variation  in  flow  angle  on 
convection  heat  transfer  ana  frictional  resistance  in  a  tube 
bank  using  a  model  scale  apparatus,  and  to  develop  constants 
which  will  be  useful  to  the  designer  of  heat  transfer  equipment. 

The  relationship  originally  developed  by  Grimiscn  [l]  *  and 
later  modified  by  the  transport  data  obtained  by  Justi  and 
Leuder L^J will  be  used  as  the  basic  convection  heat  transfer 
relationship  where  trie  tube  banks  are  inclined  at  ninety  degrees 
to  the  fluid  flow.   A  factor  F   will  be  developed  which,  when 
included  in  the  modified  version [3] of  this    1  tion  will  give 
it  tne  form: 

¥*   .  0.292  Wp  V-(   .I-"!  ,)  V3 

The  factor  F'   will  be  obtained  from  <     n  curves,  and  wnen  so 

to 

used  will  give  a  modified  value  of  Nusselt  Number  hfd,   for  the 

k 
tubes  inclined  at  an  angle  9  to  fluid  flow. 

Concurrently,  with  the  investigation  of  tne  effect  of 
variation  in  flow  angle  on  convecLion  heat  transfer  tne  effect  on 
frictional  resistance  as  expressed  by  the  Fanning  type  friction 
factor  will  be  determined. 
'Numbers  in  brackets  designate  references  listed  at  the  end  of  text 
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The  availability  of  electronic  computers  makes  possible 
more  accurate  reduction  and  analysis  of  experimental  data 
than  has  been  possible  before.   In  addition,  they  permit  a 
greater  number  of  calculations  to  be  carriea  out  in  a  short 
time.   Thus  a  systematic  series  of  geometries  will  be  tested 
by  varying  the  longitudinal  tube  pitch  and  tube  arrangement. 
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THEORY 

Effect  on  Convection  Heat  Transfer 

The  transfer  of  heat  by  a  fluid  passing  through  a  tube 
bank  is  by  the  mechanism  of  forced  convection  in  turbulent  flow. 
The  quantity  of  heat  transferred  through  tne  stagnant  fluid  film 
or  thermal  boundary  layer  contiguous  to  the  tube  is  a  function 
of  the  following  quantities: 

a)  a  characteristic  length,  taken  in  this  development  as 
the  outside  diameter  of  the  tubes,  D  . 

b)  the  coefficient  of  convection  heat  transfer,  h, ,  between 
the  fluid  and  the  tube  surface. 

c)  the  thermal  conductivity,  k,  of  the  fluid. 

d)  the  specific  heat,  c  ,  of  the  fluid  at  constant  pressure. 

e)  the  dynamic  viscosity ,>-«,  of  the  fluid. 

f)  the  mass  flow  of  the  fluid,  G. 

The  convection  heat  transfer  coefficient  may  be  expressed 

\ 

as  a  function  of  the  other  quantities: 
^   h  «  f  (Dt,  K,  cp,yw  ,  G) 

The  choice  of  a  characteristic  length  and  of  a  temperature 
at  which  to  evaluate  those  temperatures  dependent  physical 
properties  is  open  to  some  discussion.   When  fluid  flows  across 
a  single  tuoe  the  outside  diameter  of  the  tube  is  the  most 
logical  choice,  however  in  a  tUDe  bank  the  flow  is  effected  by 


17 
the  clearance  between  tubes  and  other  definitions  of  character- 
istic length  have  been  suggested. 

To  preserve  consistency  with  Pierson,   Grimison  -*  ana  Huge  ^ 
and  because  other  definitions  provide  no  proven  advantage  the 
outside  diameter  of  the  tube  will  be  used.   Similarly  various 
definitions  of  the  temperature  at  v^hich  to  evaluate  physical 
propet_ties  huve  been  suggested. 

That  which  is  consistent  with  the  work  mentioned  above 
and  which  is  most  logical  for  forced  convection  heat  transfer 
is  the  temperature  of  the  fluid  film.   As  there  is  a  variation 
in  fluid  film  temperature  through  its  thickness  the  film 
temperature  is  defined  as  the  arithmetic  mean  between  the  bulk 
temperature  and  the  tube  surface  temperature.   Likewise  the 
bulk  temperature  varies  through  the  tuue  bank  and  is  defined 
as  the  arithmetic  mean  of  the  inlet  and  exhaust  temperatures. 

The  exact  mechanism  of  forced  convection  heat  transfer  is 
not  fully  understood  so  that  while  analytical  methods  can 
sometimes  be  used,  semi-empirical  or  empirical  relationships 
are  more  reliable.   The  development  of  empirical  relationships 
is  assisted  by  tiie  derivation  of  dimensionless  groups  of  the 
physical  qualities  discussed  above. 

These  groups  may  be  derived  by  the  methods  of  dimensional 
analysis.   The  applicable  dimensions  are: 

L  for  length 

T  for  time 

M  for  mass 

H  for  heat  energy 

0  for  temperature 


As  heat  is  equivalent  to  mechanical  work  it  can  be 

2   —2 
replaced  by  < force)  x  (resistance)  or  H  *  ML   T   .   The 

dimensions  of  the  remainder  oi  the  physical  properties  may  be 

expressed  as: 


D.  =  L 

t 


h 


MT 


-31-1 


k  =  MLT^3^"1 

c   .  lV2^1 

P 

/U     =  ML*"1?"1 

-2  -1 

G  =  ML   T  x 


If  the  Buckingham  Pi  Theorem  is  applied,  the  following  matrix 
is  developed: 


M 
L 
T 


D 

0 

1 
0 
0 


h. 

1 

0 
-3 
-1 


K 

1 

1 

-3 

-1 


P 

0 

2 

-2 

-1 


JU. 

1 
-1 
-1 

0 


1 

-2 

-1 
0 


From  this  matrix  three  dimensionless  relations  whicn  nave 
proved  useful  in  heat  transfer  work  are  derived.   These  are: 

=  Nusselt  number, 


hfdt 


=  mt""3^""1  l 


-3^.-1 


MT   0 


D 


which  may  physically  be  interpreted  as  the  ratio  of  the 
characteristic  length  to  a  stationary  fluid  layer  conducting 
heat  at  the  same  rate  under  the  same  temperature  uiiference; 
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jDt   =  ML"2T""1  L      =  Reynolds  number, 

^        -1  -1 

ML   T 

which  may  be  interpreted  as  the  ratio  of  inertia  force  to 
viscous  force  and: 

-j~-    =  L2T"2j"1  .  r:L""1T~1     =  Prandtl  Number 

which  is  the  ratio  of  two  constants  of  molecular  transportation, 
The  relationship  is  usually  written  in  terms  of  power  functions 
as : 

hfDt     =  A  /GDta>\   /  C  M\h 

where  the  coefficient  A  and  the  exponents  a  and  b  are  constants, 

By  correlation  of  empirical  data  where  the  fluia  flow  is 
normal  to  the  tube  bank  the  above  equation  has  been  shown  to  be 

¥*     =°-292FAFD     /  ^t\°-60     /^/*  V3 

where  F   and  F   are  empirical  correlations  known  as  arrangement 
factor  and  depth  factor  respectively.   The  arrangement  factor 
is  a  function  of  tube  arramjenent ,  staggered  or  in-line, 
longitudinal  and  transverse  spacing  ana  Reynolds  Number.   The 
depth  factor  is  a  function  of  the  number  of  rows  in  the  tube 
bank.   The  depth  factor  determined  by  Grimison  becomes  unity 
for  ten  or  more  rows. 
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fe] 
Colburn      has    shown    that    for   fluid    flow    parallel    to    the 

tube    surface: 

hfDH      =    0.023    /    GDH   \  °'80    /C     AC 

~k-  (—J         (-f- 

Where  D,,  is  the  hydraulic  diameter.   In  oraer  to  be  able  to 
make  a  direct  comparison  with  flow  normal  to  the  tubes  the 
above  expression  may  be  written  in  terms  of  tube  diameter. 


hfDt 
k 

=    0.023 

/GDtN 

0.80 

)    ( 

'    DH 

0.20 
)             ( 

1/3 


For  cases  where  the  fluid  flow  is  neither  normal  nor 
parallel  to  the  tube  surface,  but  inclined  at  some  angle  theta 
to  it,  the  rate  of  heat  transfer  as  indicated  by  the  value  of 
Nusselt  Number  should  lie  somewhere  between  that  for  normal 
flow  or  cross  flow,  and  parallel  flow.   It  should  be  possible 
to  express  the  effect  of  the  inclined  flow  on  forced  convection 
heat  transfer  in  one  of  two  ways  either  as: 

1)   a  modification  to  the  exponent  of  Reynolds  Number,  i.e.; 
h.D,      n  OQO  _  _    rn    0.60  f(0) 

_f_t       ■  °-292  fafd/  GDt  \  /   %/u\ 


where  f  (©)  represents  a  function  of  the  angle  of  inclination,  or 
2)   as  an  additional  multiplying  factor: 

hfDt   =  0.292  FAFD  f  (©)/  GDt^0-60  f   Cp/M^   1/3 
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An  analytic  prediction  of  the  effect  of  flow  angle  on 
convection  heat  transfer  wuulci  be  at  best  uncertain.   As  the 
earlier  work  with  which  this  work  is  to  be  consistent  is  an 
empirical  correlation,  no  theoretical  prediction  will  be 
attempted.   Instead,  data  derived  from  the  test  apparatus 
will  be  correlated  witn  tne  empirical  equation  to  produce  a 
modifying  factor  which  is  a  function  of  the  angle  of  inclination. 

Effect  on  Frictional  Resistance 

The  frictional  resistance  of  a  fluid  flowing  across  a 
tube  bank  is  customarily  treated  by  an  equation  .f  the  general 
Fanning  type  wherein  the  pressure  drop  is  treated  as  Joeing 
caused  by  flow  through  a  series  of  orifices  formed  by  successive 
major  restrictions  in  the  flow  path.   By  methods  of  dimensional 
analysis  similar  to  those  described  above  for  convection  heat 
transfer  a  dimensionless  friction  factor,  f ,  may  be  derived 
which  is  a  function  of  Reynolds  Number  only: 
f  =  C  /  CD 


Where  the  coefficient  C  and  tne  exponent  X  are  constants, 
This  equation  may  be  modified  to  permit  the  use  of  convenient 
quantities  without  reducing  its  soundness: 

f  =  10.84  (10)8/:/DP 
NG2 
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where: 

10,84  =  a  constant  which  incluues  the  effect  of  the 

acceleration  of  gravity  in  ft.  per  hr.  per  hr . 

and  a  head  conversion  factor  from  pounds  per 

square  foot  to  inches  of  water 

=  air  density 

DP    =  pressure  drop  across  the  tube  bank  in  inches  of 

water 

N     =  number  of  major  identical  restrictions. 

/in  the  case  of  a  tube  bonk  N  is  equivalent  to  ; 

I  / 

i  the  number  of  rows  of  tubes 

G     =  mass  f«Low 

The  question  of  what  temperature  should  be  used  to  determine 
the  properties  of  the  fluia  again  arises.   The  temperature  has 
a  significant  effect  en  gas  density  and  viscosity  which  directly 
effect  the  friction  factor  and  trie  corresponding  Reynolds  Number. 
Utilization  of  the  mean  bulk  temperature  in  the  bank  has  been 
shown  to  be  unsatisfactory  by  previous  investigators.   Empirical 
developments  led  to  the  establishment  of  the  convention  that 
the  physical  properties  be  evaluated  at  a  temperature  equal 
to  the  tube  surface  temperature  minus  eight  tenths  the  log  mean 
temperature  difference  for  staggered  arrangements  or  equal  to 
tube  surface  temperature  minus  nine  tenths  the  log  mean 
temperature  difference  for  in-line  arrangements.   The  minus 
sign  applies  where  heat  is  transferred  from  heated  tubes  to 
a  cooler  fluid. 
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As  is  true  with  the  convection  heat  transfer  tne 
prediction  of  variation  of  friction  factor  with  flow  angle 
does  not  lend  itself  to  analytic  prediction.   Accordingly  no 
attempt  at  analytic  prediction  will  be  made. 

The  choice  of  a  purely  empirical  approach  while  not  as 
intellectually  stimulating  as  a  theoretical  development  has 
the  distinct  advantage  of  producing  results  which  are 
reproduceable  and  can  be  used  with  some  confidence  in  the  design 
of  heat  transfer  equipment. 
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EXPERIMENTAL  TECHNIQUE 

The  technique  followed  in  the  investigation  is  that  of 
forcing  air  through  a  bank  of  heated  tubes  which  are  arranged 
in  various  geometries  and  inclined  at  various  angles  to  the 
air  flow  in  order  to  measure  changes  in  parameters  indicative 
of  the  quantities  of  heat  transferred  ana  the  frictional 
resistance  in  the  bank.   Data  is  obtained  for  three  geometries, 
all  with  a  pitch  normal  to  the  air  flow  of  1,50  tube  aiameters. 
Pitches  in  the  direction  of  flow  are  1.50  and  2.00  tube  diameters 
and  the  tubes  are  arranged  in  staggered  and  in-line  configura- 
tions.  For  each  geometry,  measurements  are  taken  with  the 
tube  bank  at  angles  of  90,  60,  45  and  30  degrees  to  the  air 
.flow.   The  angle  is  varied  by  inclining  the  tube  bank  in  a 
test  section  of  constant  cross  section  and  length. 

The  tubes  are  internally  heated  by  steam,  the  inlet  pressure 
being  controlled  so  that  the  tube  surface  temperature  is  very 
near  212  degrees  Farenheit.   Ambient  air  is  heated  by  being 
forced  through  the  tube  bank  by  a  steam  turbine  driven  blower. 
Variations  in  the  mass  flow  of  air  are  obtained  by  varying  the 
speed  of  the  supply  blower.   The  mass  flow  of  air  is  measured 
by  means  of  a  square  edged  orifice  and  appropriate  manometers. 
Temperature  sensors  in  the  air  duct  upstream  and  downstream 
from  the  tube  bank  are  used  to  determine  inlet  and  outlet 
temperatures.   Tube  surface  temperatures  are  determined  by 
sensors  located  on  the  tube  surfaces  tnroughout  the  bank, 
generally  on  the  upstream  sector  of  the  tube.   Tube  surface 
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temperatures  used  in  the  computations  are  tnen  the  arithmetic 
mean  of  the  measured  tube  temperatures.   The  principle  parameter 
for  determining  the  effect  of  frictional  resistance  is  the 
pressure  drop  across  the  tube  bank  measured  by  manometers. 
Transport  data  is  evaluated  at  appropriate  temperatures  with  a 
humidity  correction  obtained  by  standard  psychometric  techniques 
applied  to  the  specific  heat. 

In  addition  to  the  aata  obtained  from  the  three  geometries 
tested,  the  data  determined  by  previous  investigators  has  also 
been  included.   All  this  aata  has  been  placed  through  the 
computation  process  and  analyzed  for  correlation  of  both  the 
previous  work  and  the  extended  investigation. 
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PROCEDURE  FOR  REDUCING  LXPERIMLKTAL  DATA 
The  purpose  of  the  reduction  of  experimental  data  is  the 
calculation  of  Nusselt  Numbers,  Reynolds  Numbers  and  friction 
factor.   This  process  may  be  divided  into  four  parts: 

1)  calculation  of  weignt  rate  of  flow  through  the  apparatus 

2)  calculation  of  Nusselt  Number  and  the  corresponding 
Reynolds  Number 

3)  calculation  of  friction  factor  and  the.1  corresponding 
Reynolds  Number 

4)  calculation  of  Nusselt  Numbers  by  the  modified  Grimison 
equation  ana  the  Colour n  equation. 

This  section  develops  in  detail  the  calculation  procedure. 

1)   Weight  rate  of  flow: 

The  weight  rate  of  flow  calculation  is  perhaps  the  most 
complex  of  those  performed.   The  equation  is: 


Vv'AO  =  A2  (AK)  Y  .  /  2g  P     (P±    -    P 2) 

where: 

WAO  =  weight  rate  of  flow 
AK  =  flow  coefficient 
Y  =  expansion  factor 
g  =  gravitational  constant 
P  =  air  density 
(P-j-P^J   =  pressure  differential  across  the  orifice 
The  method  used  to  obtain  weight  rate  of  i low  from  orifice 
dimensions  and  pressure  differential  is  a  standard  engineering 
technique  and  hus  been  developed  in  detail  by  both  Bond  ana 
Tallin  J  and  Kiss  and  Smith;   it  will  not  be  repeated  here. 
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By  combining  constants  and  applying  proper  unit  conversion 
factors  the  equation  reduces  to: 

WAO  =  o. 02415  (IDO)2  (AK)  Y  /uHF 

where: 

IDO  =  orifice  meter  diameter  in  inches 

DHF  =  pressure  drop  across  the  orifice  meter  in  inches  of  oil 

The  expansion  factor  Y  is  a  function  of  diameter  ratio, 

ratio  of  differential  pressure  to  absolute  static  pressure  and 

the  ratio  of  specific  heats.   As  shown  in  Appendix  D  the  factor 

Y  ca"n  be  calculated  for  each  aata  point  by  applying  an  equation 

of  the  form 

Y  =  1  -  C   (DHF) 

p   -  13.596PB 

The  flow  coefficient,  AK,  is  a  function  of  orifice  diameter, 

aiameter  ratio,  ana  pipe  Reynolds  Number.   The  aerivation  of 

the  expressions  for  flow  coefficient  in  Appendix  D  gives 

nations  of  the  form: 

AK  =  C2  +  C3 


/  ANRE 

However  Reynolds  Number  in  the  pipe  is  given  by  the 

expression: 

ANRE  =  Y  (IDO)    (AK)  C  /DHF 
M 

where : 

ANRE  =  Pipe  Reynolds  Number 

V  =  Expansion  Factor 

IDO  =  Orifice  Diameter 
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jj    =   Dynamic  Viscosity  for  air,  the  expression  for  which 
was  developed  in  Appendix  D  and  is   * 
^u       =    (.16  TI  +  112.6)  xl(T7 
AK  =  flow  coefficient 
C  =  a  constant 
DHF  =  differential  pressure  across  the  thin  plate  orifice 
meter. 

It  can  be  seen  that  flow  coefficient  and  Reynolds  Number 
are  mutually  dependent  and  an  iterative  process  must  be  used 
to  determine  satisfactory  values  for  these  two  quantities. 
This  can  be  done  by  first  calculating  intermediate  quantities 
composed  of  those  factors  other  than  AK  which  appear  in  the 
Reynolds  Number  expression  specifically: 


V  = 


2.239 
0.16TI  +  112 


X  10 


and : 

Y  *   1  -  C   (DHF) 
p   +  13.596PB 

A  trial  value  of  AK  depending  upon  orifice  diameter  may  be 
assumed  and  a  corresponding  value  of  Reynolds  Number  may  be 
calculated : 

ANRE  m    V  x  Y  x  AK 

This  value  of  Reynolds  Number  may  be  used  to  calculate  a 
second  value  of  flow  coefficient  which  in  turn  may  be  used  to 
compute  a  new  Reynolds  Number.   This  process  can  continue 
until  the  comparison  of  the  new  value  of  AK  with  tne  previously 
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calculated  value  shows  that  the  increase  in  accuracy  docs  not 
warrant  repetition.   This  final  value  of  flow  coefficient 
can  then  be  used  to  cdculate  weight  rote  of  flow  of  air 
through  the  apparatus. 

WAO  =  (86.9  78)  (IDO)2  (AK)  (X) 

When  the  weight  rate  of  flow  has  been  determined  certain 
other  values  may  be  determined.   These  include:   the  average 
temperature  in  the  tube  bank  in  uegrees  Rankine: 

TAVG  =  TI_+_TE    +  459#69 

and  the  specific  heat  of  wet  air  by  the  method  detailed  in 
Appendix  D,  the  log  mean  temperature  difference: 

XLMTD  =    T£  -  TI 


LOG  {   IT-TI 
e 


TT  -  TI 


ana  film  temperature: 

TF  =  /TI  +  TL N 


.0 


.  TT 


2.0 

When  these  quantities  have  been  computed  everything 
necessary  for  the  second  major  calculation  step,  the  computation 
of  Nusselt  Number  has  been  completea.   The  computation  of 
Nusselt  Number  may  then  be  undertaken. 

2)   Nusselt  Number: 

Nusselt  Number  may  be  defined  as 

Nusselt  Number  =      Q „ 

S  (XLMTD)      IT" 
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where: 

Q  =  total  heat  transferred 
Dt  =  tube  diameter 
S  =  total  heating  surface 
XLMTD  =  log  mean  temperature  difference 
k  =  thermal  conductivity  of  air 
The  relation  for  Nusselt  Number  may  then  be  expressed  as: 

ANNU  =  WAO  (CPW)  (TE  -  TI ) 

LTT(XLMTD)     (0.000024TF    +    0.01344) 

where : 

Q  =  WAO  (CPu)  (TE-TI) 
Jc  =  0.000024TF  +  0.01344 
Dt/S  =  effective  length  L  x  TT 

by  utilizing  the  interior  dimension  of  the  box,  minimum  possible 

tube  length,  and  combining  constants  the  equation  reduces  to: 

ANNU  =  ,-;A0  (CPW)  ( TE  -  TI ) 

XLMTD  (0.01158TF  +  5.4829). 

For  flow  in  which  the  tubes  are  inclined  the  surface  area  is 
increased.   The  surface  area  at  any  angle  of  inclination  is 
equal  to  the  surface  area  at  ninety  degrees  divided  Dy  the 
sine  of  the  angle  of  inclination,  thus  for  angles  other  than 
ninety  degrees  the  expression  for  Nusselt  Number  becomes: 

BNNU  =  WAO  (CPw)  (TE  -  TI )  SINE  Q     =  ANNU  ( SIN£  ^ 
XLMTD  (0.01158TF  +  6.4829) 
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This  is  the  proper  value  of  Nusselt  Number  for  inclined  flow 
and  is  the  one  result  of  the  second  part  of  the  calculation. 
It  is  now  necessary  to  calculate  the  Reynolds  Number  through 
the  test  section.   Reynolds  Number  may  be  expressed  as: 

BNRE  =  IVAO  (Dt)  , 
KM 
where: 

A  =  free  flow  area  through  the  test  section,  a  projected 

area  which  remains  constant  though  the  angle  is  varied 

s*  =  dynamic  viscosity  of  air  =  (.16TF  x  112.6)  x  10 

The  equation  reduces  to: 
BNRE  =  iVAO  x  102 


(0.046256TF  +  32.552)  . 

with  these  two  results,  Nusselt  Number  and  Reynolds  Number  the 
second  part  of  the  computation  is  complete. 

3)   Friction  factor 

For  the  frictional  resistance  calculation  some  quantities 
which  have  not  been  previously  calculated  are  required. 

E   new  definition  of  film  temperature  is  necessary.   This 
definition  is  a  function  of  tube  temperature  and  log  mean 
temperature  difference  and  may  be  expressed  in  the  form: 

TFF  =  TT  -  C5  (XLMTD) 
where  Cj.  is  either  0.800  or  0.900  depending  upon  whether  the 
tube  geometry  is  staggered  or  in-line.   A  Reynolds  Number  for 
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use  with  the  friction  factor  may  be  calculoted  using  the 
value  of  film  temperature  just  computed: 
CNRE  =  WAO  x  102 


(0.046256  TFF  +  32.552) . 
The  friction  factor  is  expressed  as: 
FRIFAC  =  10.84  (RHO)  DP  10_8 
N/WAO 

The  density  of  air  is  calculated  using  the  equation: 

P 
-^ —  =  RT  in  the  form 

-RHO 
RHO  =  (P,  +  13.596PB)  5.1815  / 

53.349   TFF  =  459.69 

When  RHO  has  been  obtained  the  friction  factor  may  be  calculated 

FRIFAC  =  1.084  (RHO)   (DP)  1Q8 

WAO2 
0.16732 

4)   Nusselt  Number  by  the  Grimison  Equation 

When  the  friction  factor  has  been  obtained,  the 
computation  continues  to  calculate  the  Nusselt  Number  by  the 
modified  Grimison  equation. 

The  Grimison  equation  is  expressed  as: 

TNNU  =  0.292  (ARGFAC)  F   (BNRE)  0-6°  /CPW/f\1/3 

As  the  test  section  is  ten  rows  deep  the  depth  factor, 
F  ,  is  unity  and  can  be  eliminated  from  the  expression.   The 
quantities  Reynolds  Number,  BNRE,  ana  specific  heat,  CPW  are 
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available  from  previous  work,  ana  analytic  expressions  in 
terms  of  temperature  exist  for  dynamic  viscosity  ana  thermal 
conductivity.   It  has  been  shewn  in  Appenaix  D  that  the 
arrangement  factor  may  be  expressed  as  a  function  of  Reynolds 
Number  for  specific  geometries.   With  proper  substitutions  the 
expression  for  Nusselt  Number  becomes: 


TNNU  =  0.292  ARGFAC  (BNRL)*5 


_  "j 


II  ..  (0.16TF  +112.6)  (1(T7) 


.C00024TF  +  C. 01344 


•  -J  3  3 


..hen  this  value  of  Nusselt  Number  has  been  determined,  a 
factor  FTHLTA  is  defined  by  forming  the  quotient  of  the 
experimental  Nusselt  Number  BNNU,  and  the  Grimison  Nusselt 
Number,  TNNU: 

FTHLTA  =  BNNU 
TNNU  . 

When  FTHLTA  has  been  computed,  the  computation  of  Nusselt 
Number  by  the  Colburn  expression  is  undertaken.   The  similarity 
in  the  Colburn  expression  and  that  for  the  Grimison  equation 
leads  to  simplifying  the  calculation  by  calculating  the  Col£>urn 
Nusselt  Number  using  the  previously  calculated  Grimison 
Nusselt  Number: 

CNNU  =  .02  3  (BNKE)0-80(DIAFAC)  /CP..><\1/3 

where  DIAFaC  is  a  conversion  factor  tne  tur-e  diameter  to 
hyaraulic  diameter. 

) °  * 6 ° ( aRGFAC )  ( CPft^  / J 

K 
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Dividing  CNNU  by  TNNU  ana  solving  for  CNNU: 

CNNU  =  .0  7877  (B.'.RC)  °  *2°  ( DIAP'AC  )  FNNU. 

ARGFAC 

This  form  is  readily  digestible  by  the  computer.   When  CNNU 
has  been  computed  the  calculation  of  one  Qota  point  is 
complete. 

When  the  calculations  have  Oeen  completed  for  all  data  points 
the  processing  of  experimental  data  is  complete  and  analysis 
can  commence. 
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ANALYSIS  OF  CONVECTION  HEAT  TRANSFER  RESULTS 

Figures  1  thru  4  are  composite  plots  showing  the 
experimental  variation  of  Nusselt  Number  with  Reynolus  Number 
for  each  of  the  four  geometries  which  have  been  tested.   The 
variation  in  slope  of  the  plots  as  the  angle  of  inclination 
changes  indicates  th^t  the  inclination  factor  F   is  not  a 
simple  multiplying  constant  but  is  an  exponential  function  of 
Reynolds  Number. 

Fft  is  defined  as  the  ratio  of  the  experimentally  determined 
Nusselt  Number  to  the  Theoretical  Nusselt  Number,  calculated 
by  the  modified  Grimison  equation,  at  a  particular  value  of 
Reynolds  Number.   Thus  for  every  experimental  point  a  value 
of  F_  as  Reynolus  Number  may  be  plotted. 

Rather  than  plot  these  values  by  hand,  the  data  processing 
facilities  at  David  Taylor  Model  Basin  were  utilized  to  enhance 
the  speed  and  accuracy  of  the  process.   An  existing  program 
was  modified  and  used  to  provide  the  needed  output.   This 

program  is  described  in  Davia  Taylor  Model  Basin  Report  No. 

r     . 
2037 L7j  and  in  Appendix  G. 

The  program  takes  as  input  the  values  F   and  the  corres- 
ponding Reynolds  Number  plus  identifying  codes  which  serve  to 
individualize  each  of  the  points.   The  output  is  the  equation 
of  three  lines  of  regression  which  are  produced  by  means  of 
the  method  of  least  squares.   In  addition,  there  is  a  graphic 
output  which  plots  the  three  lines  of  regression  which  are 
determined.   Samples  of  these  output  forms  are  included  in 
Appendix  G.   A  complete  summary  of  the  regression  analysis 
coefficients  is  shown  on  page  42. 
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MPOSI  TE    PLOT 
VARIATION    OF   NUSSELT     DUMBER 
WITH    REYNOLDS     NUMBER 
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COMPOSITE     PLOT 
VARIATION    OF    NUSSELT    NUMBER 
WITH    REYNOLDS     NUMBER 
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REYNOLDS      NUMBER       X      10 


FIG.  3 
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Of  the  three  regression  equations  fcr  each  geometry,  one 
represents  consideration  of  all  the  points,  one  is  determined 
by  all  the  points  within  two  standard  errors,  ana  the  third 
represents  the  line  determined  by  points  within  one  standard 
error,  where  a  standard  error  for  a  set  of  N  points  is  given 
by  the  relation: 


S.E. 


£      V 

'  where  Y  avg.  is   )    Li 

The  choice  of  which  equation  to  utilize  was  made  through 
inspection  of  the  SC4020  plots  of  the  regression  equations. 
In  almost  all  cases,  the  values  of  F   at  low  Reynolds  Numbers 
were  outside  the  two  standard  error  criteria,  and,  since  the 
first  three  points  are  among  those  used  to  make  a  first  deter- 
mination of  exponent,  these  points  influence  tne  phase  of  the 
regression  line  unduly.   Further,  the  data  for  values  of  Reynolds 
Number  less  than  2000  is  suspected  to  be  influenced  by  laminar 
flow.   Thus  only  that  data  which  fits  within  the  two  standard 
error  criteria  was  allowed  to  determine  the  regression  equation 
and  the  second  set  of  coefficients  were  used  for  further  analysis. 

For  design  work,  the  equations  as  shown  in  Table  I  may  be 
used  if  the  designer  is  working  at  a  specific  angle  of  inclination 
and  geometry.   Entering  with  Reynolds  Numoer  as  an  argument, 
a  value  of  F   will  result  which  can  be  used  as  an  additional 
factor  in  the  Modified  Grimison  equation. 
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For  greater  convenience,  two  design  plots  have  been 
prepared.   Figures  9  una  10  in  trie  cnclusion  section,  give 
the  factor  F   as  it  varies  with  angle  at  three  constant  values 
of  Reynolds  Number.   There  are  12  curves,  three  for  each  of 
the  four  tube  arrangement  geometries  analyzed. 

Entering  with  the  angle  of  inclination  and  the  geometry, 
the  designer  may  read  a  value  of  F   for  each  of  the  two  values 
of  Reynolds  Number  which  bracket  the  argument  he  requires. 
By  linear  interpolation,  F   can  be  particularized  to  the 
desired  Reynolds  Number. 
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TABLE  I 

Summary  of  Inclination  Factor  Regression  Coefficients 
For  Points  '/.ithin  Two  Standard  Errors 


Fe  -  A  *  B  (nre)M 


A 

B 

M 

ST«1.5DT 

45 

0.5087726 

1.056159Q0)""3 

0.6510998 

SL-1.5DT 

60 

0.6856272 

4.110530(10)~4 

0.9889 

STG 

75 

-0.6487573 

0.6318375 

0.1069 

90 

0.8435844 

1.62334Q0)"4 

0.7754198 

ST=1.5DT 

30 

0.6208538 

2.209353Q0)"5 

0.9889 

SL=1.5DT 

45 

0.7091857 

3.031199(10)"5 

0.9889 

IN-LINE 

60 

0.7860375 

4.310665(10)~5 

0.9889 

90 

1.155682 

-4.277594(10)~6 

1.00999 

ST=1.5DT 

30 

7.290824(10)" 

~2  4.&53093(10)~2 

0.2980699 

SL=2.0DT 

45 

0.7260181 

3.172470(10)"5 

0.9889 

STG 

60 

0.7262158 

1.741723C10)"4 

0.8012598 

90 

0.7653439 

7.028655(10)"3 

0.4712899 

ST-1.5DT 

30 

0.4970187 

3.12976K10)"5 

0.9889 

SL=2. 0DT 

45 

0.6867613 

2.62199K10)"5 

0.9889 

IN-LINE 

60 

0.8011853 

3.337236(10)"5 

0.9889 

90 

0.8751028 

2.385977(10)"5 

0.9889 
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ANALYSIS  OF  FRICTIONAL  RESISTANCE 

Figures  51  through  66,  included  in  Appendix  H,  are  plots 
of  friction  factor  against  Reynolus  Number  for  individual 
geometries  and  inclinations.   Several  points  pertinent  to 
subsequent  analysis  are  observable  from  these  curves.   The 
non-linear  behavior  of  the  full  logarithmic  plot  at  low  Reynolds 
Numbers  is  apparent  in  these  plots.   This  behavior  is  a 
characteristic  of  small  tube  apparatuses  that  is  consistent  with 
the  observations  of  previous  investigators  and  is  probably 
caused  by  laminar  flow  in  the  tube  bank.   In  order  to  keep  the 
analysis  as  general  as  possible  data  points  with  Reynolds 
Numbers  less  than  2500  have  been  omitted.   Additionally  it  can 
be  seen  from  Figure  55  that  the  data  for  the  thirty  degree 
angle,  5^=3   =  1.5DT,  In-Line,  is  completely  inconsistent  with 
all  the  remaining  data  and  is  therefore  omitted  from  the 
analysis.   The  inconsistenty  is  undoubtedly  due  to  faulty 
instrumentation  or  operator  error  during  the  collection  of 
this  data. 

The  friction  factor  can  be  expressed  analytically  as  a 
function  of  Reynolds  Number  by  fitting  the  data  to  an  equation 

of  the  form : 

fypo 
FRIFAC  =  CONST  (CNRE)    u 

using  the  method  of  least  squares. 

A  computer  program,  described  in  Appendix  H,  has  been 
utilized  to  compute  these  parameters  which  describe  the  friction 
factor  as  a  function  of  Reynolds  Number  for  a  particular  geometry 
and  angle  of  inclination.   The  equations  which  include  these 

I 

parameters    are    summarized    in   Table    II. 


TABLE    II 
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F    =    CONST    (CNRE) 


EXPO 


\ 

e 

CONST 

ST=SL=1.5DT 

90 

1.931187 

STAGGERED 

75 

0.445318 

60 

0.045465 

45 

0.65543 

ST-SL-1.5DT 

90 

0.261858 

In-LINE 

60 

0.037337 

45 

0.100758 

ST-1.5DT 

90 

0.431623 

SL=2.0DT 

60 

0.187797 

STAGGERED 

45 

0.070171 

30 

0.040784 

ST=1.5DT 

90 

0.297399 

SL=2.0DT 

60 

0.074656 

IN-LINE 

45 

0.017001 

30 

0.033258 

• 


EXPO 
•0.457575 

■0.045756 

■0.196590 

•0.043053 

■0.081586 
•0.119514 
•0.024007 

•0.086865 
•0.028274 
■0.089713 
•0.040784 

-0.055798 
•0.05331 
0.210122 
0.060968 
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Figures  5  thru  8  are  composite  plots  showing  the  variation 
in  friction  factor  with  Reynolds  number  at  aifferent  flow 
angles  for  particular  geometries.   The  plots  are  based  on  the 
equations  listed  above.   From  these  plots  it  may  be  seen  that 
the  set  for  on  ST  =  SL  =  1.5DT,  STAGGERED,  (Fig.  5)  is  not 
consistent  with  the  remaining  data.   The  data  for  this  geometry 
was  taken  by  previous  investigators  who  experienced  some 
difficulty  in  correlating  it.   This  previous  data  will  be 
omitted  from  further  analysis. 

From  the  equations  representing  the  variation  of  friction 
factor  with  Reynolds  number,  it  can  be  seen  that  tnere  is  no 
apparent  pattern  in  the  Variation  of  either  the  coeificient  or 
exponent  with  angle  of  inclination  for  particular  geometries. 

Friction  factors  for  Reynolds  Numbers  of  three  thousand 
through  eignt  thousand,  in  increments  of  one  thousand,  were 
computed  as  part  of  the  analysis.   Program  Values  of  Reynolds 
Number  =  3000,  5000  ana  8000  are  tabulated  below  to  indicate 
the  variation  with  flow  angle  for  particular  geometries. 

REYNOLDS  NUMBER  =  3000 

Friction  Factor 
GEOMETRY  30      45      60      75      90 

ST  =  SL  =  1'5Dt'  STG         ~    0.0464  0.2194  0.3087   0.1140 
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REYNOLDS  NUMBER  =  3000  continued 

Friction  Factor 
GEOMETRY  30      45      60      75      90 

Sfp=Sr=2.0D  ,  I.L.  -    0.0831   0.0972    -    0.1363 

IE       1 

ST=1.5DT,  SL=2.0DT,STG   0.0648  0.1312   0.1498    -    0.2153 
ST=1.5DT,  Sr=2.0DT,  I.L. 0.0542  0.0914   0.1144    -    0.1902 

REYNOLDS  NUMBER  =  5000 

ST*S  =1.5DT,STG  -    0.0454  0.2426  0.3016  0.1095 

ST=SL=1.5DT,I.L.  -    0.0821  0.1033    -  0.1307 

ST=1.5DT,  SL=2.0DT  STG   0.0668  0.1258  0.1476    -  0.2060 

ST=1.5DT,  S  =2.0DTI.L.   0.0599  0.1018  0.1176    -  0.1850 

REYNOLDS  NUMBER  =  8000 

ST=S  =1.5DT  STG             -    0.0445  0.2661  0.2952  0.1054 

ST=S__1.5DT  I.L.            -    0.0812  0.1093  -  0.1258 

ST=1.5DT,  SL=2.0DT  STG   0.0686  0.1206  0.1457  -  0.1977 

ST=1.5DT,  SL=2.0DT  I.L.  0.0575  0.1124  0.1206  -  0.1801 

From  these  values,  a  relationship  can  be  developed  wnich 

expresses  friction  factor  at  any  angle  as  a  function  of  tne 

friction  factor  at  ninety  degrees  by  simply  taking  ratios  of 

the  values  and  plotting  them:  i.e. 

F©  =  fQ 
f9© 
where  fQ  =  the  friction  factor  at  any  angle, 

f"9©  =  tne  friction  factor  at  ninety  degrees  as  expressed  by 

pypf) 

f,  =  CONST  (CNRE)      CONST  and  EXPO  being 
particularized  to  the  geometry. 


30 

45 

60 

90 

- 

0.6097 

0.7131 

1.000 

0.3010 

0.6093 

0.6958 

1.000 

0.3500 

0.4805 

0.6015 

1.000 

51 
The  values  so  derived  are  tabulated  below: 

n    =  3000  F© 

Kb 

GEOMETRY 

ST=SL=  1.5DT  I.L. 
ST=1.5DT  SL=2.0DT  STG 
ST=1.5DT  SL=2.0DT  I.L. 

NRE  =  5000 

ST=SL=  1.5DT  I.L.             -  0.6282  0.7904  1.000 

ST=1.5DT  SL=2.0DT  STG  0.3243  0.6107  0.7165  1.000 

ST=1.5DT  SL=2.0DT  I.L.  0.3388  0.5502  0.6356  1.C00 

NRE  =  8000 

ST=SL=  1.5DT  I.L.  -    0.6455    0.8688    1.000 

ST=1.5DT  SL=2.0DT  STG     0.3470   0.6100    0.7370    1.000 

S^l.50^  ST  =  2.0D„  I.L.  0.3182   0.6241    0.6696    1.000 
1       1   L  T 


Values  of  tne  ratio  F9  are  plotted  and  are  included  in  the 
Conclusion  section  as  Figures  11  through  13.   By  entering  the 
applicable  curve  with  values  of  flow  angle  the  ratio  of 
friction  factor  at  the  angle  €    to  that  at  ninety  degrees  may 
be  obtained.   For  values  of  Reynolds  Numbers  other  than  those 
tabulated,  linear  interpolation  may  be  used. 
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CONCLUSION 
Convection  Heat  Transfer 

Figures  9  and  10  represent  the  plotted  results  of  the 
equations  developed  by  the  regression  analysis  of  individual 
inclination  factors.   These  plots  show  the  variation  of 
inclination  factor  with  angle  for  in-line  and  staggered 
geometries.   Six  curves  are  plotted  in  each  case  representing 
constant  values  for  Reynolds  Number  of  2000,  5000  and  8000  and 
particular  values  of  transverse  and  longitudinal  tube  pitch. 
The  curves  can  be  used  directly  as  a  design  aid  by  applying 
the  value  of  inclination  factor  obtained  as  a  multiplier  to 
the  modified  Grimison  equation.   Values  of  the  inclination 
factor  for  Reynolds  Numbers  other  than  those  plotted  may  be 
obtained  by  linear  interpolation.   The  modified  Grimison  equation 
now  becomes: 

It  should  be  noted  that  the  variation  of  inclination  factor 
with  angle  of  inclination  is  by  no  means  linear  so  that  there 
are  regions  for  both  staggered  and  in-line  geometries  in  which 
relatively  large  changes  in  angle  of  inclination  cause  small 
changes  in  the  rate  of  heat  transfer.   The  reverse  situation 
also  can  be  obtained.   The  choice  of  which  region  is  optimum 
must  be  left  to  the  user  who  is  familiar  with  the  requirements 
and  limitations  of  his  particular  job. 
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» 

FRICTIQNAL  RESISTANCE 

Figures  11  through  13  represent  the  plotted  results 
of  the  analysis  of  frictional  resistance  for  each  of  the 
three  geometries  for  which  reliable  aata  was  obtained,   Each 
figure  contains  three  curves,  one  each  for  values  of  Reynolds 
Numbers  of  3000,  5000  and  8000.   Additionally  on  each  figure 
there  appears  an  equation  in  the  form: 

*-Hp»)(*,«)x 

This  is  the  equation  for  the  friction  factor  in  each  case. 

Two  uses  of  the  curve  are  suggested.   The  first  is  to 
calculate  an  absolute  value  of  friction  factor.   This  method 
should  be  used  with  caution  for  while  the  method  of  obtaining, 
reducing  and  analyzing  frictional  resistance  data  is  believed 
to  be  more  accurate  that  that  used  by  previous  investigators 
it  does  not  agree  closely  with  their  work.   The  second  use 
is  to  estimate  qualitatively  the  effect  of  varying  the  flow 
angle  on  the  frictional  resistance  in  a  particular  piece  of 
heat  transfer  equipment. 


56 


LQQ 

IN-Lf 

(V6 

£n 

f 

■  0.26c 

L 

r ' 

Uj 
u 

(SO 

/ 

2 

to 

'1 

40 

/ 

/ 
/ 

ft 

-J 

.pa 

/ 
/ 

N/ii  *( 

)00 

F3 

o 

n 

, 

/ 

(.  • 

ft 

i.no 

ft 

.so 

o 
<3 

£0 

o 

AO 

Nxtr  5 

D00 

-.J 

.?o 

o 
$ 

o 

■ 

i.nn 

- 

■QQ 

>^ 

- 

j 

.6,(: 

AC, 

.?n 

**£■* 

000 

c 

I 

) 

3 

0 

— — — — i 
6 

0 

9 

0 

ANGL 

E    OF 

IN(  L 

INA  T 

1  ON 

F/0.    // 


57 




i.nn 

STJSG 

GERci 
.S.--  /  5 

5 

do 

RD 

f  3 

QA321 

—  \J. 

z  N 

L/U  U  J 

o 

<3 

AO 

I/O 

to 

111 

PO 

>V=  S 

000 

■ 

ct 

-J 
<3 

o 

o 

1.0Q 

fin 

of 

o 

.60 

4  0 

5 

o 

PO 

/ 

N*/  5 

3  OP 

5 

> 

0 

Lf)0 

m 

- 

.60 

4  0 

*fc  =  3 

Q00 

.P0 

/ 

0 

1 

( 

) 

3 

0 

i 

o 

9 

0 

A/\ 

"irLE 

Of    It 

tCUN 

4  /"  /  0  A 

F/  G.  /<? 


58 


i  no 

1  N-L 

s  = 

1  NE 
IS  0T 

an 

y 

sL  = 

2.0  C, 

■ 

-Of 

.60 

f  =  0 

.297  f 

AO 

%r  <?' 

100 

20 

Uj 

0 

00 

inn 

Uj 

,«n 

o 

fiO 

.40 

NrfSQ 

00 

PC) 

hi. 

o 

0 

1 

<3 

i  on 

10 

o 
o 

$ 

.60 

40 

■  ^i 

Nnr=30 

00 

_.. 

.P0 

■Tit      ■* 

0 

( 

1 

3 

0 

6 

D 

9 

0 

ANGL 

E   or- 

1  N  C  i 

JNA1 

ION 

* 

F 

1  G.  13 

59 

RECOMMENDATIONS  FOR  FURTHER  INVESTIGATIONS 

The  apparatus,  test  procedure  and  method  of  reducing 
and  analyzing  data  in  this  thesis  represent  a  suitable  method 
of  empirically  determining  the  effect  of  flow  angle  on 
convection  heat  transfer  and  frictional  resistance  in  a  tube 
bank.   The  results  presented  are  limited  in  that  they  include 
only  four  possible  tube  geometries. 

It  is  recommended  that  the  series  be  extendea  to  incluae 
all  values  of  transverse  and  longitudinal  tube  pitches  likely 
to  be  found  in  modern  heat  exchangers.   Such  work  would  involve 
manufacture  of  new  tube  banks  and  sealing  strips.   It  is 
further  recommended  that  data  oe  taken  at  an  angle  of  seventy- 
five  degrees  both  for  other  geometries  which  may  be  tested  and 
for  those  for  which  this  angle  has  not  been  previously  used. 
The  additional  point  would  be  o±  great  assistance  in  firmly 
locating  the  upper  end  of  the  inclination  factor  curves. 

If  a  blower  of  larger  capacity  should  become  available 
the  test  data  for  all  geometries  should  be  extended  to  higher 
Reynolds  numbers.   Prospective  investigators  are  advisee  that 
the  limit  of  the  currently  installed  pressure  instrumentation 
has  been  reached  by  the  present  air  supply  system.   They 
are  also  aovised  that  the  manufacturing  processes  are  long 
and  tedious  and  that  the  condition  of  the  laboratory  steam 
plant  should  be  carefully  investigated  before  undertaking  to 
use  it. 

The  reduction  and  analysis  of  test  data  could  be  facilitated 
by  combining  the  three  Computer  programs  used  in  this  work  into 
one  program  to  reduce  and  analyze  data  in  one  computation. 
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APPENDIX  A 

Apparatus 

At  the  onset  of  trie  project,  it  was  felt  that  the  apparatus 
as  designed  by  Bond  ana  Tallin  and  later  used  by  Kiss  and  Smith 
woulu  not  be  adaquate  for  the  task  without  considerable 
modification.   The  apparatus  as  constructed  ana  available  had 
been  usea  for  one  basic  geometry  by  Kiss  ana  Smith  and  had 
an  inherent  disadvantage  in  that  it  was  difficult  to  make  any 
alteration  in  tube  geometry. 

The  design  of  a  new  test  section  was  therefore  undertaken. 
The  principles  which  guidea  the  design  are: 

1.  The  section  must  be  of  the  same  cross  sectional 
dimensions  as  the  original  apparatus  as  verified  by 
Bond  ana  Wall in. 

2.  The  section  must  have  the  capability  of  rapid  T_uj~e 
geometry  changes. 

3.  All  manufacturing  operations  necessary  for  completion 
of  the  design  must  be  feasible  with  equipment  avail- 
able in  the  Webb  shop  facility. 

4.  The  design  must  utilize  the  same  tube  banks  that  were 
prepared  for  the  previous  investigators. 

A  design  which  satisfied  the  four  stated  requirements  was 
developed.   The  configuration  is  generally  the  same  as  that 
used  in  the  earlier  experiments.   The  steam  supply  and  drain 
line  were  used  without  moaif ication,  as  was  the  air  supply 
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system  including  the  orifice  meter.   The  entire  test  section 
was  rebuilt  in  accordance  with  the  plan  shown  as  Fig.  15  in 
this  appendix. 

The  test  section  consists  of  two  rectangular  wooden  boxes, 
closed  at  one  end,  which  serve  as  the  entrance  and  exhaust 
portions  of  the  test  unit.   Connecting  these  two  end  pieces  to 
form  a  continuous  square  duct,  is  a  center  section  of  aluminum 
alloy.   This  center  section  is  designed  with  stationary  side 
plates  and  movable  Lop  and  bottom  plates.   These  features 
enables  the  tube  bank  tc  be  adjusted  to  any  angle  from  vertical 
through  60  degrees  from  the  vertical.   This  was  the  practical 
extent  of  adjustment  in  that  the  tube  banks  were  not  long 
enough  between  headers  to  permit  another  range.   This  is, 
however,  considered  to  be  a  sufficiently  wide  range,  as  it 
permits  angles  of  flow  incidence  from  30  degrees  to  90  degrees, 
at  any  practical  increment.   This  design  parameter  sized  the 
length  of  the  aluminum  center  section  at  30  inches. 

The  end  sections  are  made  up  of  3/4  inch  exterior  grade 
marine  plywood.   The  joints  are  made  up  with  glue  and  joined 
with  #8  x  lh    "  wood  screws.   The  sides  of  the  section  are 
routed  3/16"  by  3"  to  receive  the  sides  of  the  aluminum  center 
section  such  that  the  center  plate  will  fit  flush.   The 
wooden  sides  are  also  routed  at  the  top  and  bottom  edge,  3/16" 
wide  by  V  deep  to  receive  the  edges  of  the  sliding  top  plates. 
The  top  and  bottom  of  the  wooden  intake  ana  exhaust  sections 
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are  routed  to  a  depth  of  3/16  inch  to  effect  a  flush  fit  for 
the  sliding  plates  of  the  aluminum  center  section.   They  are 
also  routed  across  the  width  and  to  a  length  which  will 
accommodate  the  sliding  plates  through  the  range  of  inclination 
of  the  tube  bank.   This  range  is  shown  in  Figure  17,  the 
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drawing  of  the  section.   The  intake  section  as  a  6  inch  90 
elbow  secured  to  the  end  for  delivery  of  air  through  a  6M 
diameter  hole.   The  exhaust  section  has  a  similar  stove  pipe 
elbow  made  up  to  the  side  of  the  section.   The  heated  air 
exhausts  through  a  6  inch  diameter  hole  in  the  box  then  through 
the  stove  pipe  exhaust  line  to  the  atmosphere. 

Both  intake  and  exhaust  sections  were  sanded  smooth  and 
finished  witn  two  coats  of  shellac  prior  to  assembly.   After 
assembly,  all  edges  were  caulked  with  a  silicone  rubber  sealant 
and  the  stove  pipe  connections  were  maoe  up  with  asbestos  cloth 
gasketing  using  Permatex  compound  for  a  leakproof  joint. 

The  center  portion  of  the  test  section  is  fabricated  from 
2024-T3  aluminum  alloy.   This  alloy  was  chosen  for  its  free 
machining  characteristics,  since  several  machining  operations 
were  required  in  manufacturing  the  numerous  sealing  strips 
between  the  tube  banks.   It  was  later  learned  that  the  temper 
of  this  plate  was  not  correct  in  that  when  the  plates  were  cut 
to  size  and  the  edges  ana  sliding  grooves  milled,  the  plates 
warped  out  of  plane  due  to  locked  in  roiling  stresses.   After 
several  unsuccessful  attempts  to  straighten  the  plates  by 
different  methous,  the  decision  was  maue  to  anneal  the  plates 
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the  plates  and  thus  stress  relieve  them.   Several  industrial 
contacts  were  made,  ana  through  the  courtesy  and  kinaness 
of  the  Sylvania-Corning  plant  in  Hicksville,  N.Y. ,  the  plates 
were  annealed.   Provided  thus  witn  plane  surfaces,  construction 
of  the  center  section  may  proceed.   One  further  note  could  be 
added  for  the  benefit  of  those  who  may  be  curious,  a  more 
proper  alloy  choice  for  this  application  is  2024-T53  this 
temper  oeing  stress  relieved  by  stretching.   The  siae  plates 
were  made  up  with  all  edges  milled  and  fastened  to  the  wooden 
end  pieces  by  flat  heat  wood  screws,  courtersunk  to  effect 
a  plane  surface.   The  top  movable  plates  were  maae  up  having 
all  edges  milled  and  two  3/16"  wide  by  1/16  inch  deep  groove 
milled  in  the  sides  to  effect  a  seal  with  trie  siae  plates  and 
to  provide  a  channel  for  guidance  of  the  moving  plates.   Each 
top  and  bottom  plate  (four  total)  has  a  projection  along  each 
edge  which  slides  in  a  channel. 

The  channel  was  mace  from  an  aluminum  storm  window  shape 
and  proved  to  have  insufficient  strength  to  hold  the  weight  of 
the  tube  bank  and  secure  the  plates  in  position  as  well.   To 
correct  this,  heavy  steel  yokes  were  made  from  angle  and 
threaded  to  close  the  section  and  secure  the  sealing  strips  in 
place.   All  but  one  of  the  channels  were  removed,  the  one 
remaining  only  to  serve  as  a  positioning  device.   Manufacture 
of  the  sealing  strips  represented  the  largest  expenditure  of 
time  for  any  one  operation.   For  each  of  the  12  combinations 
of  angle  and  geometry,  22  sealing  strips  were  required,  11 
each  for  top  and  bottom.   Each  sealing  strip  required  edge 
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milling,  groove  milling,  and  22  holes  drilled  or  milled  to 
close  tolerances  for  completion.   In  all,  5680  holes  were 
machined,  and  a  total  oi  264  strips  were  prepared.   The  sealing 
strips  proved  very  satisfactory  in  maintaining  the  tube  geometry 
while  providing  excellent  sealing  characteristics.   The  only 
modification  required  was  to  provide  a  small  hole  in  the  end 
of  each  strip  to  permit  the  separate  pieces  to  be  wired  to 
hold  them  in  place  while  changing  tuioe  geometry.   An  isometric 
sketch  showing  the  arrangement  of  the  center  test  section 
and  the  sealing  strips  is  provided  as  Figure  15.   Drawings 
showing  the  sealing  strip  as  a  part  of  the  test  section  and  in 
detail  appear  as  Figures  17  and  18.   A  photograph  of  the 
sealing  strips  in  detail  is  available  as  Figure  22. 

The  section  is  supported  from  a  laboratory  table  by  two 
wooden  frames  and  a  framework  of  steel  angle  which  was  modified 
only  to  facilitate  the  changing  of  the  geometry.   This 
modification  consisted  of  welding  an  angle  on  each  corner  of  the 
metal  structure  to  provide  support  for  the  tube  bank  headers 
when  changing  the  sealing  strips. 

The  tube  bank  is  that  used  by  previous  investigators. 
It  is  made  up  of  10  rows  of  k    inch  O.D.  copper  tubes,  22  tubes 
per  row,  with  each  row  of  tubes  being  fitted  with  an  inlet 
and  exhaust  header  of  1%   inch  copper  pipe.   A  detail  drawing 
of  an  individual  tube  bank  is  shown  in  Figure  16.   The  inlet 
headers  are  connected  to  the  steam  supply  moisture  trap  by 
3/8  inch  copper  tubing.   Many  of  these  lengths  of  tubing  had 
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had  to  be  replaced  to  permit  the  freedom  of  movement  required 
in  the  tur>e  bank.   The  exhaust  headers  were  connected  in  the 
same  manner  to  a  drain  header  of  1%   inch  copper  pipe.   An 
additional  problem  was  that  these  lines  kinked  when  movea  about 
due  to  strain  hardening  and  several  were  replaced  for  this 
reason.   For  some  geometry  changes,  the  exnaust  and  supply 
lines  were  completely  disconnected  to  avoid  this  problem. 

Two  small  holes  are  provided  in  one  side  of  the  center 
test  section  for  thermistor  leads.   Four  holes  are  provided 
for  pressure  taps  in  the  test  section.   Two  are  in  the  inlet 
section,  positioned  in  the  middle  of  the  vertical  height  and 
4  inches  from  the  center  section.   Two  taps  are  provided  in 
the  exhaust  section  on  a  similar  vertical  position,  4  inches 
from  the  test  section.   Thermistor  probes  are  provided  in  the 
6  inch  inlet  piping  just  before  the  entrance  elbow  and  in  the 
6  inch  outlet  section  to  avoid  radiation  effects  from  the 
tube  bank. 

Previous  investigators  had  used  various  methods  of  sealing 
the  banks;   References  4  and  Q.      The  method  used  in  this  series 
of  tests,  though  long  and  painstaking  in  preparation,  proved 
worth  the  effort  while  testing.   There  was  no  need  to  use  any 
sealing  compound  once  the  tubes  were  in  place  and  the  top 
and  bottom  moving  plates  made  tight.   The  tolerance  to  which 
the  sealing  strip  holes  were  prepared  was  0.008  inch  and  there 
was  little  or  no  leakage.   Further,  a  geometry  change  could 
be  made  in  about  two  hours  on  the  average  and  there  was  no  glue, 


69 


sealer,  etc.  to  clean  from  the  tuce  to  cause  further  delay. 
In  this  regard,  the  extra  effort  was  rewarded  and  this 
arrangement  is  recommended  for  later  investigators. 

Figure  14  shows  an  overall  isometric  view  of  the  steam 
plant  and  test  apparatus. 
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Instrumentation 
Temperature 

Instrumentation  for  obtaining  temperature  readings  was  by- 
thermistors  installed  in  air  intake  and  exhaust  lines  ana  on 
tubes  to  sample  temperatures  throughout  the  tune  bank. 

Thermistors  are  solid  state  electronic  semi-ccnauctors 
which  exhibit  a  large  negative  temperature  coefficient  of 
resistance.   Their  small  physical  size  and  the  large  negative 
temperature  coefficient  of  resistance  make  thermistors 
particularly  well  suitea  for  application  in  small  scale  exper- 
imental models  as  units  can  be  located  within  a  very  limited 
space  and  great  accuracy  can  be  obtained  with  a  relatively 
simple  electrical  circuit. 

The  thermistors  used  were  Fenwal  GB35J1  bead  type  thermistors. 
These  units  are  0.043  inches  in  diameter  and  have  a  nominal 
resistance  of  5K  ohms  +  20%  at  25   C.   The  thermistors  were 
calibrated  by  the  use  of  a  constant  temperature  bath.   Boiling 
water  and  boiling  benzene  were  used  to  obtain  calibration 
beyona  the  range  of  the  bath.   A  sample  calibration  curve  is 
included  as  Figure  23  at  trie  end  of  this  section.   Thermistors 
used  for  inlet  and  outlet  temperatures  were  fitted  into 
probes  made  from  brass  tuning  and  sealed  with  silicone  rubber 
sealant.   The  probes  were  located  in  the  center  of  the  inlet 
and  exhaust  ducts  at  a  point  around  a  ninety  degree  bend  from 
the  tube  bank  so  thdt  they  would  be  screened  from  direct 
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radiation  from  the  tube  bank.   Some  difficulty  was  experienced 
in  affixing  the  thermistor  units  tightly  \.o    the  tu~e  bank,  so 
that  a  true  tube  temperature  would  be  indicated.   The  problem 
was  finally  overcome  by  taping  the  thermistors  to  the  tuDe 
with  plastic  electricians  tape.   Prior  to  placing  the 
thermistor  on  the  tube,  a  small  groove  was  ground  on  the  tube 
surface  to  receive  the  beaa.   The  tape  and  pigtail  leaas  were 
then  covered  witn  a  high  density  epoxy.   This  effectively 
held  the  thermistor  in  place  even  curing  the  curing  period. 
Thermistor  surface  sensor  assemblies  consisting  of  thermistor 
beads  bonded  to  teflon  or  silicone  impregnated  fiberglass  tape 
are  commercially  available.   The  cost  of  these  sensor  assemblies 
is  approximately  eight  tim  s  that  of  trie  thermistor  units 
so  chat  it  was  not  considered  economically  feasible  to  use  thera. 
For  similar  reasons  it  was  decided  to  undertake  calibration 
of  the  units  rather  than  to  purchase  units  which  track  a 
standard  curve  within  narrow  limits. 

The  thermistor  leads  were  connected  to  a  selector  switch 
which  switched  individual  units  into  a  Wheatstone  Bridge  circuit. 
The  resistance  necessary  to  balance  the  bridge  was  measured 
and  this  value  was  used  to  enter  the  calibration  curves  to 
determine  temperatures  in  degrees  Farenheit.   A  diagrammatic 
sketch  of  the  temperature  measuring  circuit  appears  as 
Figure  24. 

A  high  failure  rate  was  experienced  with  the  thermistor 
units.   The  hignest  percentage  of  failures  were  due  to 
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accidental  severing  of  the  small  diameter  wire  pigtails  at  the 
point  where  the  wire  entered  the  glass  beaa.   These  failures 
were  caused  simply  by  the  rough  service  to  which  the  units 
were  subjected,  particularly  those  installed  on  the  tube 
surface.   In  addition,  several  thermistors  were  found  to  be 
erratic  in  behavior,  culminating  in  a  short  circuit  within 
the  bead  after  a  short  time.   No  explanation  for  this  performance 
has  been  found.   Faulty  units  were  replaced  with  spares  until 
the  supply  was  exhausted.   At  this  time,  it  became  necessary 
to  replace  the  exhaust  temperature  thermistor  with  a  mercury- 
in-glass  thermometer. 

As  an  aside,  it  should  be  added  that  the  thermistors  were 
not  given  fair  trial  in  this  application.   The  necessity  to 
move  the  tube  bank  about  roughly  placed  undue  strain  on  these 
units.   Based  on  this  experience,  haruy  commercially  prepared 
mounts  such  as  those  mentioned  above  are  recommended. 
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Pressure 

Instrumentation  for  pressure  wus  by  manometers  containing 
oil  or  water  as  appropriate. 

The  manometer  used  to  determine  the  pressure  on  tne 
upstream  side  of  the  orifice  was  an  open  water  manometer.   The 
barometer  reading  was  converted  from  inches  of  mercury  to 
inches  of  water  and  added  to  the  manometer  reading  to  obtain 
an  absolute  pressure  upstream  of  the  orifice. 

An  oil  filled  U-tube  difierential  manometer  was  used  to 
indicate  the  pressure  drop  across  the  orifice  plate.   A  water 
filled  U-tube  differential  manometer  and  a  slant  tube  manometer 
were  used  to  measure  the  pressure  drop  across  the  tube  bank, 

A  Bourdon  tube  pressure  gauge  and  a  bellows  type  compound 
gauge  were  used  to  measure  steam  inlet  ana  drain  pressures 
respectively. 

Flow 

Instrumentation  for  weight  rate  of  flow  was  by  thin  plate 
orifice  meter  with  pressure  tops  located  at  one  pipe  diameter 
preceeding  and  one  half  pipe  diameter  following  the  inlet  face 
of  the  orifice  plates.   To  obtain/a  large  spectrum  of  Reynolds 
Numbers,  three  separate  orifices  of  2  inch,  3  inch,  and  4 
inch  diameters  were  used.   The  weight  rate  of  flow  was  obtained 
from  the  differential  head  across  the  orifice.   It  was  found 
that  it  was  possible  to  cover  the  entire  range  of  weight  rates 
of  flow  which  could  be  obtained  with  the  available  steam 
turbo-blower  with  the  3  inch  orifice. 
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APPENDIX  B 

Test  Equipment  Operating  Procedure 

Operation  of  the  test  equipment  can  be  divided  into  two 
categories:   1)   Lining  up,  lighting  off  and  operation  of  the 
steam  plant  and  auxiliary  equipment  and  2)  Operation  of  the 
test  apparatus  itself. 

The  following  procedure  is  followed  in  putting  the  steam 
plant  into  operation: 

1)  The  cooling  tower  is  checked  and  filled  as  necessary 
to  ensure  a  proper  quantity  of  cooling  water. 

2)  Cooling  tower  fans  and  the  circulating  pump  are 
started. 

3)  The  condensate  pump  is  started. 

4)  Fires  are  lighted  in  the  boiler  following  the  procedure 
given  in  the  boiler  operation  and  maintenance  manual. 

5)  When  the  steam  pressure  has  reached  about  100  psi  trie 
air  ejectors  are  lined  up  and  sufficient  nozzles  cut 
in  to  maintain  bbout  ten  inches  of  vacuum  in  the 
condenser.   The  most  satisfactory  operation  was  obtained 
when  both  second  stage  nozzles  and  no  first  stage 
nozzles  were  used. 

6)  The  turbo  blower  is  lighted  off  using  a  throttle  pressure 
of  about  5  psi.   The  blower  lubricating  oil  tank  should 
be  checked  to  ensure  that  an  adequate  supply  of 
lubricating  oil  is  flowing. 
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7)   The  steam  separator  is  drained  and  tne  tube  bank 

steam  supply  and  drain  valves  are  opened.   A  pressure 
of  ten  psi  at  the  steam  separator  was  found  to  be 
very  satisfactory. 

When  the  above  operations  have  been  completed  and  the 

conditions  in  the  apparatus  as  reflected  by  inlet,  exhaust, 

and  tube  temperatures  have  stabilized  the  procedure  enters  its 
second  phase,  the  collection  of  data. 

Test  procedure  consists  primarily  of  varying  the  quantity 
of  amounts  of  air  supplied  to  the  apparatus  and  recording  the 
data  provided  for  by  the  installed  instrumentation,  as  described 
in  Appendix  A. 

The  weight  rate  of  flow  of  air  to  the  test  section  is  varied 
by  throttling  tne  blower  supply  steam  with  fine  control  being 
accomplished  by  manipulation  of  the  discharge  Joleed  valve. 
Satisfactory  operation  is  obtained  by  starting  at  low  flow  rotes 
and  increasing  gradually  to  high  rate..   With  tnis  procedure, 
the  apparatus  stabilizes  quickly  so  that  the  interval  oetween 
runs  is  very  short.   It  should  be  noted  that  because  the  test 
section  and  the  turbo  blower  are  supplied  by  steam  from  tne 
same  line  fluctuations  in  boiler  pressure  caused  by  variations 
in  blower  demand  will  cause  fluctuations  in  the  pressure  to  the 
tube  banks.   As  the  tube  temperatures  and  exnc.ust  temperatures 
are  sensitive  to  the  supply  pressure  variations  should  be 
minimized  by  ensuring  the  tube  bank  supply  is  the  same  before 
each  run. 


The  pressure  drop  across  the  orifice  is  indicative  of 
flow  rate,  accordinqly  successive  runs  may  be  maue  by  adding 
small  increments  of  pressure  drop  by  throttling  the-  blower. 

During  testing,  it  is  advantageous  to  make  a  running 


plot  of  the  quantity  TT-TI  yDHF   versusyDHF  on  full 

TT-TE  J 

logarithmic  paper.  TT-TI  \J  DHF  may  be  considered 

TT-TE  . 

directly  proportional  to  Nusselt  Number  while ^/ DHF  may  be 

considered  directly  proportional  to  Reynolds  Number.   A  linear 

plot  on  full  logarithmic  paper  then  serves  as  a  check  on  the 

consistency  of  the  data.   Care  should  be  taken  not  to  take 

such  a  plot  as  indicative  of  the  validity  or  accuracy  of  the 

data  as  the  proportionality  is  not  exact.   A  sample  plot  appears 

as  Figure  25  at  the  end  of  this  Appendix. 

Orifice  plates  of  three  different  diameters  are  available 
for  use  with  trie  apparatus.   Utilization  of  trie  three  inch 
orifice  permits  data  to  be  taken  over  the  entire  range  of 
flow  obtainable  with  the  turbo  blower.   The  other  two  orifice 
sizes  were  not  used. 

If  an  immediate  change  of  geometry  is  anticipated  upon 
completion  of  uata  taking,  allowing  the  blower  to  run  after 
the  steam  supply  to  the  tube  bank  has  been  secured  cools  the 
metal  portions  of  the  apparatus  quickly  and  facilitates  handling. 

Securing  procedure  is  essentially  the  reverse  of  the 
lighting  off  procedure.   Upon  securing  in  cold  weather  the 
cooling  tower  should  be  drained  to  prevent  freezing  of  the  pipes. 


88 


RUNS  2601-2626 


p: 

-  : 

' 

— T— I- 

— r- 

T v  :  !    : 

r 

1 

: 

■  .•'  • 

i 

-j- 

TXi 

- 

:    ' 

*- 

-H 

... 

— -- 



i 

' 



— — 

—  - 

rrrt:  r~ 

- — 





_  — 

: 

■   :    , 

— 

r 

-1 

■?' ;  • 

'    .   ; 

, 

— 

i 

.     j. 



- 

.      •  .j 

, 

-          ■ 

. 



" 

— _ 

"' 

i 

1 
1 

4-  - 

i 

I 

LSdl 

j 

- 

1 

! 

m 



' 

i  • 

■ 

•  . 

■ 

•  - 
t~t — — ■ 

-— 

■ 

. 

!^ 

-           H 

i 

r  '    j 

. 

— pr- 

^~  q — ! 

-■        -';    . 

_1 

\] 

y 

y* 

X 

_^ — 

l 

I    - 

JfL— 



-^>i 

i    . 

■  — 

—     | 

i 

yr' 

'""""'. 

i  1 

h>' 

: 

: 

- 

- 

•L: 

■ 

\— 

i ..  . 

;- 

- 

© 

1 

' 

i . 

i  _fr 

i 

■ 

i 

-   i  -' 
i 

1 

.  .  , 

i— 

1 

• 

• 

' 

. 

1 

1 

. 

' 

""r 

■ 

" 

■ 

■ 

1.5 

2 

2.5 

3 

4            5          6 

7 

8      9 

in 

1.5 

2 

2.5 

3 

i 

3 

6 

7 

8       9     1 

/DHF 


FIG    25 


89 

APPENDIX  C 
Development  of  Equations  Used  in  Computations 

Computerization  of  the  calculation  procedure  required  the 
development  of  equations  for  a  number  of  quantities  which  are 
usually  tabulated  or  taken  from  curves  in  engineering 
applications.   Development  of  equations  which  are  applicable 
in  the  range  of  temperatures  ana  pressures  which  exist  in  the 
apparatus  reduce  the  complexity  of  the  program  and  permit  a 
great  saving  in  core  storage  within  the  machine.   The  quantities 
for  which  equations  are  developed  herein  are: 

1)  Thermal  conductivity  of  air 

2)  Dynamic  viscosity  of  air 

3)  Specific  heat  of  wet  air 

4)  Flow  coefficients  for  orifice  meters 

5)  Expansion  factors  for  orifice  meters 

6)  Grimison  Arrangement  factors 

7)  Nusselt  Number  by  Colburn  equation 

The  symbols  used  in  this  appendix  are  those  developed  for 
use  in  the  computer  program, 

1)   The  thermal  conductivity  of  air  varies  in  an  essentially 
linear  manner  in  the  range  of  temperatures  which  occur  in  the 
apparatus.   An  equation  may  be  written  in  the  form: 

k  =   /k 

t 

Where 


x  T  =    C 


k'  =  thermal  conductivity  of  air  in  Btu 

hr-ft  °F 
Ale  =  change  in  thermal  o  nauctivity  per  aegree  Farenheit. 

At 


so 


T  =  Temperature  in  degrees  Farenheit 

C  =  a  constant 
Values  of  thermal  conductivity  used  to  develop  the  equation  were 
taken  from  Table  2  of  Gas  Tables  by  Keenan  ana  Keys.  [9 J  The 
equation  for  thermal  conductivity  is: 

k  =  .000024  (T-40)  +  .0144 
or 

k  =  0.000024T  +  0.01344 
This  equation  gives  values  which  correspond  exactly  witn  the 
tabulated  values  between  40  and  240  degrees  Farenheit. 

2)   Dynamic  viscosity,  like  thermal  conductivity,  varies 
linearly  with  temperature  in  the  range  of  application  and  can  oe 

expressed  in  the  same  form: 

><  -  AM 

AT   X  T  +  c 

Where : 

/*   =  dynamic  viscosity  in  Ids 

sec-f t 

-  change  in  dynamic  viscosity  per  degree  Farenheit 

At 

T  =  Temperature  in  degrees  Farenheit 
C  =  constant 
Values  used  in  the  development  of  the  equation  were  taken  from 
Table  2  of  Gas  Tables  by  Keenan  and  Keyes.   The  equation  for 
dynamic  viscosity  is: 

yU    =  (".16  (T-40)  +  11'.    x  10" 


-7 


-7 


/<.   =    f.lST  +  112.6    x  10 
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3)  The  development  of  equations  for  specific  heat  of 

wet  air  at  constant  pressure  requires  the  selection  and  adoption 
of  applicable  equations  for  the  specific  heats  of  air  ana  water 
and  combining  them  with  proper  weighting. 

The  equations  selected  were  taken  from  Table  I  of 
Thermodynamics  by  F aires  [lOj. 

For  dry  air: 

CPD  =  0.219  +  0.342T  x  10"4  -  .293T  x  10"8 
and  for  water: 

CPH  -  1.102  -  33.  IT""'*  +  416 

T 

where:  T  =  average  temperature  in  the  tube  bank  in  degrees 

Rankine  the  specific  heat  of  wet  air  is  then  determined  by: 

CPW  =  SPHU  (CPH)  +  7000  CPD 
SPHU  +  7000 

where:   CPW  =  specific  heat  of  wet  air 

CPD  =  specific  heat  of  ary  air 

CPH  =  specific  heat  of  water 

SPHU  =  specific  humiuity  in  grains  per  pound 

4)  Flow  Coefficients 

For  square  edged  thin  plate  orifice  flow  meters  with  taps 
one  pipe  diameter  upstream  and  one  half  pipe  diameter  downstream 
the  following  equation  for  flow  coefficient  obtains: 

AK  =  Ko  +  b 
where : 

AK  =  flow  coefficient  corresponding  to  any  particular 
set  of  the  parameters:   pipe  diameter  (D),  ratio  of  orifice 
diameter  to  pipe  diameter  (3,  ,  and  pipe  keynolus  number  (ANRE) 


92 


Ko  =  Flow  coefficient  when  pipe  Reynolds  number  becomes 
infinite 
b  =  an  empirical  multiplier 
=  1000 


\/anre 

Ko    and   b   may   be    expressed    in    the    following    forms: 
(0.6014    -    0.01352    D~"*    +    (0.3760    +    0.07257   d"""'1 

0.00025 
D2fi2      +    .0025D    +     £4      +    1.5^16 
b    -    (0.0002         0.0011)    +    (0.0038    +    0.0004)    X 


D 


D 


\B2    +    (16.5    +    5D/516) 


For  situations  such  as  exist  in  the  test  apparatus  where 
there  are  a  limited  number  of  possible  values  of  pipe  diameter 
and  diameter  ratios  specific  values  of  Ko  and  b  may  be 
calculated  and  applied  to  the  basic  equation.   It  should  be 

■J  C. 

noted  that  although^    represents  a  number  less  than  unity 
raised  to  a  relatively  large  power  it  cannot  be  ignored  if 
acceptable  accuracy  is  to  be  obtained. 

The  values  for  Ko  and  b  for  the  various  orifice  plates 
available  are: 

Two  inch  orifice: 
Ko  =  0.59781 

b  =  .00080 
Three  inch  orifice: 
ko  =  0.61S0 

b  =  0.00133 
Four  inch  orifice: 
Ko  =  0.6735 

b  =  0.00229 
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The  expressions  for  flow  coefficients  for  each  orifice  then  become 
Two  inch  orifice: 

AK  =  0.59781  +  0.80 

AN  RE 

Three  inch  orifice: 

AK  =  0.6180  4  1.33 
v'ANRE 

Four  inch  orifice: 

AK  =  0.6735  +  2.29 

Janre 

5)   Expansion  Factor 

The  equation  for  the  expansion  factor  through  a  square 

edged  thin  plate  orifice  meter  witn  taps  one  diameter  upstream 

and  one  half  pipe  aiameter  downstream  may  be  expressed  as: 

V  =  1  -  (0.41  +  0.35/  4)  X_ 

K 

where:     •'  =  ratio  of  orifice  diameter  to  pipe  diameter 

X  -    ratio  of  differential  pressure  to  absolute  static 

inlet  pressure 

K  =  ratio  of  specific  heats. 

4 
If  a  constant  value  of  K  =  1.4  xs  assumed  and  values  of/.: 

are  calculated  for  the  two,  three  and  four  inch  orifice  meters 

the  equation  reduces  to: 

For  2"  orifice:   Y  =  1  -  .295X 

For  3"  orilice:   Y  =  1  -  .307X 

For  4"  orifice:   Y  -  1  -  .339X 

For  utilization  in  the  program  the  factor  X  must  be 
expressed  in  the  same  units  as  the  input  data.   The  numerator 
X  is  equal  to  the  differential  pressure  across  the  orifice 
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meter  (DHF)  which  is  in  inches  of  oil.   The  denominator  is 

equal  to  the  sum  of  the  barometric  pressure  (PB)  expressed 

in  inches  of  mercury  and  the  pressure  above  barometric  ahead 

of  the  orifice  meter  (P-j)  expressed  in  inches  of  water.   All 

units  will  be  converted  to  inches  of  water  so  the  expression 

for  X  becomes: 

X  =  0.82  7DHF 

P   +  13.596PB 

where:  DHF,  P  ,  and  PB  are  input  quantities  as  defined  in  the 

preceeding  paragraph. 

0.827  =  tne  specific  gravity  of  the  oil  used  in  the  manometer 

13.596  =  the  conversion  factor  from  inches  of  mercury 

to  inches  of  water. 

Inserting  the  expression  for  X  and  simplifying  the  expressions 
for  expansion  factor  as  used  in  the  calculation  are  obtained: 
For  2"  orifice: 

Y  =  1  -     0.24464  DHF 

PjL  +  13.596PB 

For  3"  orifice: 

Y  =  1  -     0.25457  DHF 

P   +  13.596PB 

For  4"  orifice: 

Y  =  1  -     .28131  DHF 

P   +  13.596PB 

6)   Grimison  Arrangement  Factors 

The  arrangement  factors  used  in  the  Grimison  Equation 
are  a  function  of  longitudinal  ana  transverse  pitches,  Reynolds 
Number  and  the  arrangement  (staggered  or  in-line)  of  the  tube 
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bank,  and  are  not  amenable  to  a  simple  analytic  expression. 
However  as  only  one  transverse  pitch,  1.5  tuoe  diameters,  and 
a  limited  number  of  combinations  of  back  pitch  and  arrangement 
were  tested  it  was  found  feasible  to  plot  arrangement  factor 
as  a  function  of  Reynolds  Number  only  for  each  geometry.   The 
curves  so  plotted  were  parabolic  in  form.   An  equation  of  the 
form  C   (BNRE)   +  C2  (BNRE)  +  C3  =  ARGFAC  could  be  written  for 
each  of  them  by  utilizing  the  LaGrange  Method.   The  equations 
so  developed  are  of  the  form: 

ARGFAC  =  C   (BNRE)2  +  0.        (BNRE)  +  C3 

where : 


GEOMETRY          C 

l    x    (10)1C 

C2  x  (10)5 

C3 

ST  =  SL  =  1.5  Dt 

IN-LINE 

.463 

-.2963 

1.063 

ST  =  S.  =  1.5  D, 
L           t 

STAGGERED 

-1.29 

.3462 

.94059 

ST  =  2.0  Dt,  SL  = 
IN-LINE 

1.5  Dt 
5.7407 

-1.774 

1.10518 

ST  =  2.0  Dt,  SL  = 
STAGGERED 

1.5  Dt 
4.0  74 

-1.5  74 

1.1698 

Nusselt  Number  by  the  Colburn  Equation 

In  order  to  establish  a  lower  limit  wnich  would  be  comparable 

to  parallel  flow  the  Colburn  equation  was  used  to  calculate 

Nusselt  Number.   The  Colburn  equation  is: 

,1/3 


hfDh   =  .02  3  /  GD 


K 


h 


M- 


where : 


K 


GD    =  Reynolds  Number  based  on  hydraulic  diameter. 

For  the  Colburn  Nusselt  Number  to  be  directly  comparable  with 
these  developed  experimentally  or  by  the  Grimison  equation  it 
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must  be  based  on  tube  diameter  rather  than  hydraulic  diameter. 
To  accomplish  this  both  sides  of  the  equation  are  multiplied  by 


the  ratio  of  tube  diameter  to  hydraulic  diameter. 

0.80 


hfDH   D 


=  0.023/v:; 


1/3 


K 


hfDfc 


D 


1! 


=  0.02  3  /G 


in  the  nomenclature  used 


CNI.'U  =  0.02  3  (BNRE) 


0.80 


n 


0.20 


D 


k  -' 


1/3 


H 


The  hyuraulic  diameter  is  defined  as  four  times  the  free  flow 
area  divided  by  the  wetted  perimeter.   For  a  staggered  tube 
arrangement : 


DH  =   4 


STSL   -  TTDt 


2T4-) 


! 


JT  D. 


DH  =   4 


STSL  "  ilDt 
TT  D, 


5 


\ 


j—1 


A 


\ 


^\:, 


For  the  in-line  arrangement: 


D 


H 


DH  ■ 


7 


L 


'T^L    "^  T 


TT  D 


1 1 


'1WL 


7TD. 


7" 


j — i 


or 
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Since  only  one  value  of  ST  and  two  values  of  S..  are  considered, 
the  hydraulic  diameter  can  be  numerically  evaluated  for  these 
values.   For  S_,  =  1.5D,  ,  S-.  =  1.5D 


DH  =  4  (1.5D^)2  -  D^TT 


DtTT 


(2TT- 


Similarly  for  St  =  1.5Dt,  Dl  =  2.0Dfc 


(i^ji) 


The  factor/D, \  '    then  can  oe   evaluated  for  a  particular 

geometry.      This  factor  has  been  given  the  name  DIAFAC  and  is 
included  in  the  Colburn  equation  to  give  it  the  form: 

CNNU  =  .023  (BNRL)0,80  (DIAFAC)  /ct>wM1/3 

*   J 
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APPENDIX  D 
Experimental  Data  Reduction 
Program  Description 

The  computer  program  described  herein  is  a  FORTRAN  IV 
program  designed  for  use  with  the  IBM  7090  data  processing  system 
The  program  calculates  Nusselt  Number,  Reynolds  Number  and 
frictional  resistance  coefficient  for  varying  flow  angles  from 
experimentally  derived  aata  in  the  Haeberle  Laboratory,  /'.ebb 
Institute.   The  constants  ana  to  some  degree  the  method  used  are 
r ticularized  to  the  apparatus  used  for  data  collection. 

The  first  step  in  the  program  is  to  calculate  the  sines  of 
the  angles  of  inclination  at  which  tests  were  conducted  and  to 
assign  to  them  floating  point  variable  names  in  order  that  they 
may  be  repeatedly  called  without  recalculation  of  the  sine 
function. 

Control  of  the  program  is  achieved  tnrougn  a  DO  loop 
utilizing  the  variable  M  which  represents  the  sequential  number 
of  a  set  of  aata  being  processed.   It  is  required  that  the 
maximum  value  of  M  called  KASES  be  supplied  as  input  each  time 
the  program  is  used.   The  value  of  M's  incremented  by  one  each 
time  a  set  of  aata  is  processed.   The  DO  statement  refers  the 
machine  to  a  READ  statement  causing  a  new  set  of  data  to  be 
read  in  and  processed.   When  the  maximum  value  of  M  (M=  KASES) 
is  reached  the  program  is  completed. 

Required  input  to  the  program  is: 

1)   Run  code  numbei ,  a  four  digit  fixed  point  variable, 
(IRCM) 
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2)  Pressure  drop  across  the  thin  plate  orifice  meter 

in  inches  of  oil  a  three-  or  four  digit  floating  point 
variable,  (DHF). 

3)  Pressure  above  atmospheric  ahead  of  the  orifice  meter 
in  inches  of  water,  a  three  or  four  digit  floating 
point  variable,  (PI), 

4)  Barometric  pressure  in  incnes  of  mercury,  a  four  digit 
floating  point  variable,  (PB). 

5)  Pressure  drop  across  the  tube  bank  in  inches  of  water, 
a  three  or  four  digit  floating  point  variable,  (DP). 

6)  Temperature  of  the  inlet  air  in  degrees  Farenheit,  a 
three  digit  floating  point  variable,  ( TI ) . 

7)  Temperature  of  the  exhaust  air  in  degrees  Farenheit, 
a  four  digit  floating  point  variable,  (TE). 

8)  Temperature  of  the  tuoe  surfaces  in  degrees  Farenheit, 
a  four  digit  floating  point  variable,  (TT). 

9)  The  diameter  of  the  orifice  plate  in  inches,  a  single 
digit  fixed  point  variable,  (IDO). 

10)  The  specific  humidity  of  the  inlet  air  in  grains  per 
pound,  a  three  digit  floating  point  variable,  (SPHU). 

11)  A  geometry  factor  indicating  whether  the  tube  geometry 
is  staggered  or  in-line,  a  single  digit  fixed  point 
variable ,  (LGF) . 

12)  An  angle  factor  indicating  the  angle  of  inclination  of 
the  tube  bank,  a  single  digit  fixed  point  variable, 
(LANGLE) . 

13)  Back  pitch  factor,  indicating  whether  the  back  pitch 

is  1.5D  or  2.0D,  single  digit  fixed  point  variable  (LAF) 
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The  program  first  calculates  the  intermediate  quantity 
V  for  use  in  the  pipe  Reynolds  Number  calculation. 

The  program  then  queries  the  value  of  orifice  diameter  to 
determine  the  path  to  be  taken  in  calculating  the  Reynolds 
Number  thru  the  orifice.   The  constants  chosen  for  calculating 
the  expansion  factor  and  flow  coefficients  are  functions  of  the 
diameter  of  the  orifice  meter.   The  procedure  followed  is 
the  same  for  all  orifice  diameters.   A  value  of  flow  coefficient 
(AK)  is  assumed,  and  a  value  of  V  is  calculated. 

The  three  quantities  which  have  been  calculated  are 
multiplied  together  to  determine  an  initial  estimate  of  Reynolds 
Number. 

When  an  estimate  of  Reynolds  Number  has  been  maae  a  better 
estimate  of  flow  coefficient  is  obtained. 

The  estimate  of  flow  coefficient  is  then  compared  to  the 
initial  value.   If  the  two  values  vary  by  one  half  of  one  per 
cent  or  less  the  calculated  value  of  BK  is  taken  as  the  value 
of  the  flow  coefficient  and  the  value  of  Reynolds  Number  is 
accepted  and  the  calculation  continues.   If  the  difference 
between  the  two  values  of  flow  coefficients  is  greater  than  one 
half  of  one  per  cent  the  calculated  value  of  flow  coefficient 
is  used  for  an  initial  value  for  recalculation  of  Reynolds 
Number  and  flow  coefficient.   This  process  continues  until  the 
one  half  of  one  per  cent  criterion  is  reached. 

When  satisfactory  values  of  Reynolds  Number  and  flow- 
coefficient  have  been  established  the  program  proceeds  by 
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squaring  the  orifice  diameter  and  converting  it  to  a  floating 
point  variable  for  use  in  the  weight  rate  of  flov;  calculation. 
Weight  rate  of  flow  is  then  calculated. 

The  average  temperature  in  the  tube  bank  is  calculated 
and  converted  to  degrees  Rankine.   With  this  average  temperature 
the  specific  heats  of  dry  air  and  water  are  calculated  and 
combined  to  determine  the  specific  heat  of  the  air  in  the 
apparatus. 

Log  mean  temperature  difference  and  film  temperature  are 
then  calculated. 

All  quantities  necessary  for  the  calculation  of  Nusselt 
Number  are  now  available  and  line  quantity  is  calculated. 

The  angle  factor  (LAi.'GLE)  is  used  to  select  the  proper  value 
of  the  sine  of  the  angle  of  inclination  previously  calculated. 
This  sine  function  is  then  multiplied  by  the  calculated  Nusselt 
Number  to  determine  the  Nusselt  Number  for  the  angled  configuration, 

The  Reynolus  Number  in  the  tube  bank  is  calculated. 

The  Program  then  queries  the  geometry  factor  to  determine 
if  the  geometry  is  staggered  or  in-line  and  branches  accordingly 
to  calculate  film  temperature  for  the  friction  factor 
computation. 

Reynolds  Number  for  use  in  the  friction  factor  calculation 
is  calculated. 

The  density  cf  the  air  in  the  apparatus  and  the  friction 
factor  are  calculated. 
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After  the  friction  factor  has  been  computed  the  program 
again  queries  the  geometry  factor  and  then  queries  tne  back 
pitch  factor  to  determine  the  geometry  so  that  tne  proper 
arrangement  factor  and  diameter  factor  for  use  in  the  Grimison 
and  Colburn  equations  can  dq    computed.   when  the  arrangement 
factor  has  been  selected  ana  computed  it  is  employed  in  the 
Grimison  equation  to  calculate  a  vdlue  of  Nusselt  Number.   The 
value  of  Nusselt  Number  calculated  by  the  Grimison  equation  is 
divided  into  that  obtained  by  use  of  the  experimental  uata. 
value  of  Nusselt  Number  is  calculated  using  the  Colburn  equation. 

All  computations  have  been  completed  so  that  an  output 
tape  may  be  written.   Output  consists  of: 

1)  Run  Code  Numoer  as  read  in  (IRCN) 

2)  Nusselt  Number  determined  from  experimental  aata  (BNNU) 

3)  Nusselt  Number  determined  from  the  Grimison  equation  (T  NU) 

4)  Nusselt  Number  determined  from  the  Colburn  equation  (CJMNU) 

5)  The  ratio  between  the  two  Nusselt  Numbers  (FTHETA) 

6)  Reynolds  Number  based  on  the  film  temperature  used  in 
the  Nusselt  Number  calculation  (BNRE) 

7)  Friction  factor  (FRIFAC) 

8)  Reynolds  Number  using  the  film  temperature  used  in  the 
friction  factor  calculation. 

Upon  completion  of  writing  the  output  data  another  set  of 
data  is  read  in  and  the  program  repeated.   With  tne  7090  system 
calculation  of  each  set  of  data  .equires  about  one  seventh  of 
one  second. 
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A  block  diagram  of  the  program  is  presented  herewith.   The 
block  diagram  is  a  graphic  representation  of  the  data  processing 
system  procedures  used  in  the  solution  of  the  problem.   The 
diagram  is  included  to  assist  the  reader  with  a  visual  record 
of  operational  and  logic  steps  in  the  program  and  the  relaLionship 
each  step  in  the  program  has  to  the  solution. 

The  symbols  used  in  the  diagram  are  the  standard  IBM  symbols 
representing  machine  functions.   A  Key  to  these  symbols  is 
included  in  the  Nomenclature  section. 

A  program  listing  also  follows  the  block  diagram. 

The  computer  facilities  used  for  these  computations  were 
those  at  the  David  Taylor  Model  Basin.   Data  was  prepared  using 
equipment  available  at  Webd  Institute. 


BASIC       DATA       COMPUTATION   

COMPUTER    PROGRAM     BLOCK     DIAGRAM 


CALCULATE 
TRIAL    RP/NDllDS 
NO.  AND  FlDio 
r.rFFF/riFAiT 


RK?QK 


com  pure' 

RATIO 
C--  AH/ BK 


COLUMN     HEADINGS 
FOR     OUTPUT 


NUMBER    OF  DATACARDS 


CALL    FOP\ 
SINE  TUNC  TKO/V 
FOR   ^0,tS°.Ln° 
\f\N.O   75°    / 


LOOP-  DO  FRO, 
I  TO  NUMBE, 
OF  DATA  CARP. 


RUN  CODE  NUMBER,  DHF  Pi,  PB ,  DR  TT,  TE.  TT 
ORIFICE  SIZE,  SPECIFIC  HUMlDlFX,  GEOMETRY 
FACTOR,   ANGLE  FACTOR,  BACK  PITCH   FACTOR 


CALCULATE 
V 


2 


/LO&IC 

Select  o/?i  f/C/ 


A5    LEFr 


AS   LEFT 


BK<AK 


AK-'&X 


DIFFERENCE     BETWEEN  ASSUMED 
AND  CALCULATED     VALUES    OF  FLOU/ 
COEFFICIENT 


COMPUTE 

RATIO 
C«  BK/ 'AS 


C  <  O.  996 


AH'  BK 


LANCLE  -5 


HUSSELT  NO 

T I  FIE  5 
SINE  JO ' 


dumm  y   5/7)  tement 
to  srtisfy  system, 
requirements 


■Q 


TOLERANCE   LIMIT  FOR  ASSUMED    AND 
COMPUTED    VALUES    OF  FLOW  COEFFICIENT 


SET  &;i Si 
Fl.i  .v  COEFFI 
AAIP     REITS, 


VI-! ED 


CALCULATE 
WEIC/4T  FLOW 
MERAG/Z  TET1, 

-,H&  Iff    Njr^r 


REI  TERF. 
C>   0.99S 


TONS   CONTINUE  UNTIL. 


CAICULSFE'  /LOG 
MEAN  TEMP  Dl, 
MEAN  TEMP 


LMGlE'--  <o.7.9 


NUSSFLT   NO 

TIME  5 
5INE  415° 


LAN£££± 


LAN -J. 


MULTIPLICflT/ON    OF /VUSSEZ.T 
NUMBER  BY  THE  J  WE  OF  THE 
AM  _-  _.-    r.,~  /  ,jCi_,  V/J7-/OK     TO 
HC COUNT  FOTr    /NCRENSSO 
TUBE  SURFACE    F)RE/1 


I  ANGLE  =  7 


NUSSELFNO. 

TIMES 
SINE  tkO° 


NUSSELT  NO. 

T/MES 
SINS'   76" 


BNNU  = ANNU 


CALCULATE 
REYNOLDS 
NUMBER 


CONTINUE 


■OT^Wf    TZ/3ES 


STAG  GERED 


CHLCULAFE 
FILM  TEMP 
TT-O.8  LMTD 


CRLCULHTE 
FIL  M     TEMP 
TT-  O.  SLMTD 


LE 


RE 


CALCULATE 

VNOLD'o    NO. 
NS/TY,  FRICTION 


t 


DUMM  Y 


CONTINUE 


DUMMY 


CONTINUE 


3 


CALCULATE 

\rRAN-5-FMEM 

FACTOR 


5FL  EC  r 

D/A  MLTTEA? 

FACTOR 


NUS 
BY 


Z*£Cl- 

u&y/ 

Si_= 

I.CTDr 

Si~ 

Z.OUr 

CALCULATE 
'IRRANOECMEA 
FAC  TOR 

T 

SAM  £ 

AC 

LEFT 

SELECF 
Dl  4  MECTBT? 
FAC  TOfR 

5AM  £ 

AS 
LE  FT 

CALCULATE 

>ELT    NUMBE 
CR/R/ISOU  FQ CATION 


R 


CALCULATE 
FTHETA 


» 


NIL 

e 


CALCULATE 
SSELT  NUMB 
COL  BURN  CO 


'f? 

JATION 


YES 


SFOP 


J 


3/        / 


i   A  ^  I 


FN 


SOUI  iTATFMI 


IF  (  M 


FCf  IRAN  PRHtiiAN'   (  0  CALCGLAT    NU  :    L  T   JUM   R.  REYNOLDS  NUN    !<   AND 
IFKICTIC.   I  ACTOR   IN  A  CONVECTICN  '  ANK  USING  I  XPERIMENTAL  DATA 

1  FCRVAT( I) 

2  FCR^AK  1H1,  3X,  32FiRUN  fV  COL/3X  •  78HCOD!  NUS 
1SFLT             nUSSFLT            MUSSEL!            FTHETA                     YNCLD           F     tICTION        REYNOLDS< 

)         ) 

3  f  ORi'  AT< 14.3F7.2.F7. 3.F6.1 .2F7.1 , ]    '  .         .  I  .      14,13) 

4  Ffil     'Al  (    'X,  r  4,  2X  ,F8.  4,  2XW8.  4,      '.         .4,2X,F      .      ,  2X,I       .      ,2X«F9«6*2X, 
1.2) 

SIN     =S  !(■•■<      .       :36  ) 

IN  i-S  I     (      .  7654  ) 

S  I N  6  =  S  I  .M  (  1  .     4  7  2) 

5  I  N  7  -  S  I   ,  (  I  .  3  ) 

A  i .      <   :  ,  1  )      K  A  S 
.  I  = 
19      !  (  M  -  i  ,  K  A  ;  t 

IF  (  f-  T  )  i         ,  I         ,17 

i       -•■  -  I  ;     (     ,  ?  ) 

^T~     ••• 
17    Kl   AC  (  5,      )  IRC.N,  DI-.F  •  P  i  , PO,DP,  I  I  ,  1      ,  T  T  ,  I  DO  •  3PHU.LGF  »L  ANGLE  ,LA<-; 

V  -  <      .23  9/  ((C.16*TI)+112.6))*<1.     L-  7  ) 
IF(  ICO-      ) 5, 6,  7 

<5     t  K  =      . 

Y-(l.  -(   .24464*DHF )/<Pl+13.596*PR> )*SGRT(DHF) 
14  AnR  -AK «7*Y 

-  ,     3 1  +  0.8:1  /SQRTC  ANRE  ) 
GC  I  ' 

6  £K=     .' 

Y  -  <  1  .     -  (  .  •  2 5457 *C  HF  )  /  (  P  1  +  1  3  .  5  S 6  * P  b  )  )  * SOR T  {  DHI    ) 
It     AKR;    =AK«V*Y 

6K  =  "  .61         +1  -33/S'JRT(  ^NRi    ) 
CC 

7  A  v  _      . 

Y-  (  1  .  .-(      .28131  »     ;-if   )/<Pl+13.596*P8) >*S  iRT(OHF) 
16     ANR         A* k V*Y 

t      .       i/SOR F( ANR      ) 

8  IF(  AK-  >9,  1      ,11 
11    C=   -;</. A  K 

CO      !"r       1- 

9  C  =  Ak/.-> 

12:     IF(      .99'.-  C  )  1     ,1      ,13 
13     AK=,  h 

IF( IDO-3) 14» 15, 16 

i :    a  k = 

=  i  DC*  ■ 
WAG-     <•   .      78*^  I  *A!    a  Y 
TAVG=(T I+TE )/2.      f 459. 69 

CPC=      »219+<      ..'4.»TAV'V1.         4)-(      .293*TA  VG**2/1  .  :  ) 

CPh-1.1     2-<  33.  l/SQRK  TAVG)  >  ♦(  ^16.     /TAVG) 
CPW-(  <SPHU*CPH)+(  70  ::     . " *CPD) ) /< SPHU+71         .      ) 
XLVTi  =(  r     -TI  )/ALOG{  (  TT-T  I  >  /  (  TT-  Tf    )) 
TF=((TE<TI)/2.C+TT)/2. 

ANNU=(WaU*CPW*<  Tt-T  I  )  )/(  XLV  TDM      .      11    >8*TF  +  6.4  82     )  ) 
IF (langl; -4 ) ?4,        .     . 
24     ENNU=ANNU*S1N3 


'/     / 


v  m  r 


FN 


IUHCE      STATFMI     i  f        -         ^  (•(  5) 


GDI         3 

25  ;  NNCJ  -  ANi\>U*S  1 i\4 
G(      1(1     3 

26  IF(LANGL  lf-7)27,  28,29 

27  FjNMJ  =  ANmU*5  IN6 

^e    £NNU=AN\U*S IN7 

GC     TfJ     3 
29    t.NNU=AN  •  U 

LNRt=WAu*<  1  •     i    2  )/(  (      .     4,       56     TF)+3     .        ■      ) 

EF(LGF)         , 2  1.22 

21  TFF=1 T-(      .8        *XLMTD ) 
GC    TO     2  3 

22  TFFrTT-l      .9...     *XLMTD ) 

23  CNRfc=(*'Ai..i*(  1  .  ))/((.     462:i6*TFF  )+32.552) 

HO=( (Pl+13.5  >6*P6)«      .  1      15  )/(  L :3  .    149*  (TFF  1-4  59.  (  9)  J 
t      I  k  AO  (  1  mC    >4  *RI  D*DP*{  1  .   »F8  )  )  /{  (  >  AG/     .  1673  )**;    ) 
IF(LGF),      ,31,32 

31  IF(LAF  )         ,33,  34 

32  IF(LAF)Z     ,3 '3, 

33  ARGFAC=(4...    74*(  L<NRfc**2)  )/(  1  .OE1      )  -  (  1  .  574*  (  BNK>|    ))/(!.         5)+1.1698 
1    PIAFAC=     .883 

bC     !■'     3  7 

34  ARGFAC=(5.74     7*  (  ^NRf**  2  )  )  /(  1  •     r"  1      )-(  1  •  7  74  »())/(  1  •     E5) +1 • 1 "51 ( 
341     DIAFAC=     .  S278 

C  0     r  f]      3  7 

35  ARGFAC=      .94     5*  + ( 3 • 4&2*BNRE ) /l .         ■-(      .129*(8NRI    **      ))/i..       - 
351     CIArAC-      .833. 

g  c    r  n    57 

36  ARGF  AC-(4.63*(H  **2  )  )  /  1  •'  fZ  1  1-  (  2  •  963*BNKE  )  /l  •  "         U,      ^6 
361     CIAI-AC-     .82 

37  TNNb=( UNPE**     .5     )*f,  .  392*ARGFAC*  (  (  <CP'**3.6*(      •  1  6*TF+  1 1  2.  6)  )  /  (  (  1  •        4 
1)*(      .  ,P4*TF+     .      1344)))**     .333) 

Fit-,     ra-      MNU/TNNU 

38  CNNU=(TN(\U*     .     7t"77*(BNRE»*      .  >     )  *D1  AF  AC  )  /AttGFAC 

.ill»(6,4)      IKCN,hNNU,TNNU,CNNU,FTHETA,i      IRF.I     JIFAC.CN 
MT  =  !"    T-l 
CCNTINU, 
oTOP 
ND 
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APPENDIX  E 
Experimental  Data 

This  Appendix  presents  a  listing  of  all  test  data.   The 
data  is  listed  sequentially  by  run  code  numbero.   No  lobels 
are  provided  on  the  data  listings,  accordingly  the  following 
information  is  provided: 

Column  Description 

1  Run  code  number  for  data  point  identification, 
four  digits,  (RNCD) . 

2  Pressure  drop  across  the  thin  plate  orifices, 
inches  of  oil,  three  or  four  digits  (DHF) 

3  Pressure  ahead  of  the  thin  plate  orifice,  inches 
of  water,  three  or  four  digits  (P-,) 

4  Barometric  pressure,  inches  of  mercury,  four 
digits  (PB) 

5  Pressure  drop  across  the  tube  bank  inches  of 
Mercury,  four  digits  (DP) 

6  Inlet  temperature,  degrees  Farenheit,  four 
digits  (TI) 

7  Exhaust  temperature,  degrees  Farenheit,  four 
digits  (TE) 

8  Tube  surface  temperature,  degrees  Farenheit, 
four  digits,  (TT) . 

9  Orifice  diameter,  inches  (IDO) 

10  Specific  humidity,  grains  per  pound  (SPHU) 

11  Geometry  indicator  'LGF) 

12  Angle  indicator  (LAN'GLE) 

13  Arrangement  indicator  (LAF) 


1401 
.402 
.403 
LA04 
L405 
[406 
L407 
L408 
1409 
L41G 
[411 
1412 
1413 
[414 
.41? 
[416 
[417 
.416 
.419 
[420 
1421 
,422 
.423 
.474 
.425 
[426 
.427 
[426 
[429 
.430 
.431 
.432 
[433 
[434 
,435 
,436 
[60] 
.602 
,603 
[604 
L605 
,606 
.607 
L608 

[609 
,61f 
[611 

612 
.613 
.614 
[615 
.616 

617 
,618 
.619 


1.4 

1. 

1.7 

1.7  . 

2. 

2.00 

2.50 

2iJ50 

2.70 

2.80 

2. 

3. 

3.90 

4  •  1  0 

4.70 

4  .  3  ( 
4. 
4.9. 

5  •TO 
6.20 
6.5< 
7.80 
8.3., 
8.30 
8.5 
9.60 

11. A: 
13.  lo 
16.00 
16.1; 
1 8  •  4  ( 
2  ).50 
24.00 
25.5. 
29.K. 
33. 
1.3 
1.60 
2.6- 
2.9< 
3.3. 
4.1 


6.00 
6 .  i> 
7.00 
7.20 

.1 

•AC 

■ 
•  6C 

.60 

.60 


8 

8 

9 

10 

11 

11 
12 
12.90 


1. 

1.40 

1.20 

1.    V 

1.20 

1.60 

1.20 

1.80 

1.60 

1.90 

2  •  Q  0 

2. 

2.40 

2.50 

3.00 

3.20 

3.00 

3.00 

3.60 

4. 

4. 

5.00 

6.20 

5.2  0 

6.40 

6.50 

7.60 

9.20 

9.80 

10.50 

12.20 

13. 

16. 

16.80 

19.20 

22. 
.50 
1.2C 
1.60 

2.40 

2. 

2.60 

3.30 

4.00 

4.2f 

4.60 

5. 

5*44 

5.4  0 

6.40 

6.6  0 

7.40 

6.60 

8.60 

8.60 


3  0.46 

29.70 

3^.46 

29.70 

29.95 

3  0.46 

30.46 

29.70 

2  9.95 

29.70 

30.46 

29.95 

29.70 

29.95 

30.46 

2  9.7: 

2  9.9  5 

30.46 

29.95 

29.70 

30.46 

29.70 

29.70 

30.46 

2  9  .  7 

30.46 

29.7' 

30.46 

29.70 

30.46 

30.46 

30.46 

30.46 

29.70 

30.46 

30.46 

30.48 

30.4  8 

30.48 

29.62 

30.48 

30.46 

30.48 

30.  4'; 

29.82 
30.48 
29.82 
29.82 
30.48 
30.48 
29.82 
29.82 
30.48 
.48 
2  9.,. 


.340 

0.120 

0.100 

0.160 

0.130 

0.120 

0.490 

0  .13  0 

0.160 

0.130 

0.160 

0.210 

0.210 

0.710 

0.800 

0.230 

C.230 

0.310 

0.360 

0.570 

0.390 

1.550 

0.390 

1.630 

0.570 

0  .  6 

0.6  70 

0.830 

0.83  0 

0.830 

1.040 

1.220 

1.300 

1.400 

1.4 

0  •  2  00 

0.400 

0.500 

1.100 

0.600 

0.900 

1.2  00 

1.400 

1»9 

1.700 

2.1 

2.4  00 

1.700 

2.500 

2.C 

3.2  00 

2.900 

4.< 

3.600 


80.0 
70.0 
78.0 
79.0 
72.0 
72.0 
78.0 
72.0 
72.0 
79.0 
72.0 
72.0 
79.0 
73.0 
73.0 
72.0 

?3.0 
7A.r, 

74.0 
79.0 
74.0 
79.0 
78.0 
75.0 

so.o 

78.0 
8  0  •  0 
80.0 

82.0 
S2.0 
63.0 
64.0 
6  5.0 
84.0 
87.0 
68.0 
74.0 
75.0 
76.0 

81.0 

76.0 
77.0 
76.0 

78.0 

82.n 

78.0 

81.0 

82.0 

78.0 

30.  0 
82.0 
83.0 
61.0 
83.0 
63.0 


176.0 
166.0 
174.  r 
176.0 
165.0 
163.0 
174.0 
162.0 
163.0 
167.0 
163.0 
164.0 
176.0 
163.0 
163.0 
162.-' 
162.- 
162.0 
165.' 
166.0 
160.0 
167.0 
15  5.0 
15  7.0 
157.0 
160.0 
162.0 
161.0 
160.0 
159.0 
156.0 
159.0 
156.0 
156.0 
156.0 
155.0 
171.0 
164.0 
165.0 
166.0 
165.0 
164.0 
161.0 
160.0 
162.0 
156.0 
161.0 

161. 0 

15  8.0 

155. C 

160. 
160. "> 
156.-: 
154.0 

159.0 


712.0 

712.0 
212.0 
212.0 
212.0 
212.0 
2  1. 2  .  0 
212.0 
212.0 
212.0 
217.0 
212.  0 
212.0 
212.0 
717.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.r, 
212.0 
212.0 
712.0 
712.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
212.0 
217.0 

212.0 

212.0 

212.0 

212.0 

212.0 

212.0 

212.0 
212.0 

7 1 ; . 

212.0 
212.0 

212.0 
212.0 
212.0 
212.0 

212.0 


3 
4 
3 

3 

3 

3 

4 
3 

3 
3 
3 
3 

3 
3 

4 

3 
3 
3 
3 
3 
4 
3 
4 
3 
3 
3 
3 
3 

3 
3 
3 
3 

3 

3 
3 
3 
3 
3 
3 
3 
3 

3 
3 
3 

3 
3 
3 
3 
3 
3 
3 
3 


77.0 

BT.O 

77.0 
94.0 
70.0 
77.0 
77.0 
87.0 
70.0 
94.0 
77.0 
70.0 
94.0 
70.0 
77.0 

87.0 
70. 0 
77.0 
70.0 
94.0 
77.0 
94.0 
87.0 
77.0 
87.0 
77.0 
04.0 
77.0 
94.0 
77.0 
77.0 
77.0 
77.0 
94.0 
77.0 
77.0 
62.0 
62.0 
62.0 
46.0 
62.0 
62.0 
62.0 
62.0 
46.0 
62.0 
46.0 
46.0 
62.0 
6  7.0 

46.0 
46.0 

62.0 
62.0 
46.0 


0 
0 
0 

0 
a 

0 

0 

o 
o 
o 

0 
0 

o 
0 

0 
0 

0 

0 
0 

o 

0 
0 

0 

0 

d 

o 

0 
0 

o 
0 

o 

0 

o 
0 
0 

o 

r. 
0 
o 
0 

o 

o 

0 

o 

0 

o 
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APPENDIX  F 
Reduced  Experimental  Data 

This  Appendix  consists  of  a  listing  of  the  results  of 
the  data  reduction.   The  column  headings  are  self  explanatory, 
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APPENDIX  G 
Convection  Heat  Transfer  Regression  Analysis 
A  regression  analysis  is  applied  to  the  results  of  the 
data  reducing  computer  program,  the  purpose  of  which  is  to 
fit  a  curve  of  the  generalized  type,  Y  =  A  +  BX   through 
the  reduced  data.   The  independent  variable  X  is  the  exper- 
imental value  of  the  Reynolds  Number  through  the  test  section, 
defined  throughout  this  text  as  BNRE.   Tne  independent 
variable  Y  is  tne  ratio  of  the  Nusselt  Number  as  measured 
experimentally  to  the  Nusselt  Number  computed  from  the 
modified  Grimison  equation,  defined  by  FTHETA.   The  results 
of  this  regression  analysis  will  be  the  constants  A  and  B, 
and  the  exponent  M  resulting  in  an  equation  relating  FTHETA 
to  Reynolds  Number  BNRE,  within  the  range  of  the  experimental 
data,,   From  this  equation,  particularized  for  the  angle  of 
inclination  and  the  tube  geometry,  the  angle  factor  FTHETA 
may  be  computed,  or  taken  from  design-oriented  curves  to  apply 
as  an  additional  factor  to  the  modified  Grimison  Equation. 

The  analysis  was  carried  out  using  a  prepared  program 
available  at  the  Applied  Mathematics  Laboratory  of  the  David 
Taylor  Model  Basin.   The  program  is  written  in  Bell  FORTRAN  II 
for  the  IBM  7090  computer.   Running  time  is  2.5  minutes.   The 
program  is  documented  in  DTMB  Report  20  3  7. 

A  brief  description  of  the  program  is  included  in  this 
Appendix. 
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A  maximum  of  35  aata  points  may  be  entered  as  input  for 

any  one  case.   The  program  will  fit  a  function  of  the  form 

M 

Y  =  A+B  X   to  the  data  points  entered  by  the  method  of  least 

squares.   The  line  of  regression  will  be  calculated  and  plotted 
three  times  using  three  sets  of  data  points,  the  first  being 
the  entire  original  set  of  data  points,  the  second  being 
the  set  of  points  within  two  standard  errors,  ana  trie  third 
using  those  points  which  fall  within  one  standard  error.   The 
standard  error  is  defined  as: 


SE  = 


Z  {  Vl  v1_"  ' 


'--(     N 

where:  Y  average  is  the  average  ordinate  of  all  the  points. 

Yi  is  the  general  ordinate  when  i  is  indexed  from  1  to  N, 
N  is  the  number  of  points. 

The  method  used  by  the  program  is  the  method  of  least 
squares.   An  approximation  to  the  exponent  M  is  made  by  using 
the  first  three  and  the  last  three  points,  averaging  the 
ordinates  and  abcissas  of  these  points  to  derive  an  average 
point,  X  average.   Four  points,  (two  greater  than  and  two  less 
than)  closest  to  X  average  are  taken.   Frcm  these  four  points, 
a  new  average  point  is  derived.   The  slopes  of  the  upper  and 
lower  portion  of  the  first  approximate  line  of  regression  are 
then  determined.   Diviuing  the  lower  slope  by  the  upper  slope, 
an  expression  exists  from  which  the  exponent  M  may  be  approximated. 
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x,   K 


*«•>  YL 


XS  =  Xi+X2+X3 


Y   =  Y  +Y  +Y 
S     12   3 


L     n   n-1  n-2 


YT  =  Y  +Y   _+Y   0 
L     n   n-1   n-2 


X  avg. 


■  VXL 


The  four  values  of  X  closest  to  X  avg.  are  located 


X-,  =  X.  .  +X.  .  +  X4  +  X.  +  1 
C     i-2    l-l     i     l 


Yn    =  Y.  0  +  Y.  -  +Y.  +Y. 

C     i-2     l-l    i    l+l 


in  which  X.  =  X  avg. 


Then  the  exponent  M  may  be  determined  in 
the  following  manner: 

M 

Y  =  A  +  BXJ 


4*   =  BMX*1-1 

dx 
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(i)  vc-ys  (*c   ♦  xs 


Two  slopes  may  be  written  from  Fig.  A 

M-l 


XC  ~  xs 


\ 


where  Y_  -  Y   is  the  slope  at  X„  +  Y^ 


XC  ~  XS 


YT  -  Yr  /  XT  +  Xrx^ 


where  Y   -  Yc  is  the  slope  at  X   +  X„ 
XL  "  XC  ~ 


Dividing  equation  (2)  by  equation  (1). 

(Y   "  Yc)  (Xc  -  Xs)  /XL  ♦  X^1 


(X_  -  X  )  (Y_  -  Yc)  VX-,  +  X 

L     q  C      b  \  C 


Now  let   (Y   -  Yc)  (Xc  -  Xs)     Q 
(XL  *  V  (YC  -JYS) 


Let   XT  +  X_ 

L     C  =  E 

xc  +  xs 

D  =  E  <M-X) 


LnD  =  (M-l)  In  E 
Then  M   —    In  D 


c 


+  1 


In  E 
which  provides  an  approximation  to  M. 

Tne  initial  approximation  of  the  exponent  is  testing  for 

accuracy  by  considering  seven  values  of  the  expon.nt,  M  being 

the  center  value  and  the  others  1/4  M  from  this  central  value  in 
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increasing  ana  decreasing  increments. 

For  each  value  of  the  exponent,  the  values  of  A,  ana  B 
are  computed,  ana  the  sums  of  the  squares  of  the  derivations 
are  minimized.   The  exponent,  M,  is  then  chosen  from  the  set 
of  A,  3,  and  M  that  has  yielded  the  minimum  sum  of  the  deriva- 
tions squared. 

M 
The  normal  equations  for  the  function  Y  =  A  +  BX   with 

A,  B,  and  M  as  the  regression  parameters. 


nA    +    B    £  X .  -    1 

*-     l 

A      X.           +    B  Z.    X. 

l  l 


.  Y.    =    0 

2M       -2    Y.X.M 

1    1 


=    0 


A  I.     in    X.     XM      +    Blx.2M 
li  l 


M 


In    X.     -  > Y.X.       In      X.     =    0 
l  li  i 


These  equations  are  solved  for  A  ana  B  by  Cramer's  Rule. 
The  final  test  for  the  proper  exponent  is  maae  by  examining 
the  coefficients  of  the  three  normal  equations,  i.e.,  three 
equations  in  two  unknowns.   Since  A  and  B  are  unique,  two  of 
the  three  equations  must  be  linearly  dependent  ana  the 
determinant  of  the  coefficients  will  be  identicaily  equal  to 
zero. 
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2M 


Y.X. 
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i  X 


M 


iX 


2M 
X.    In  X. 
i        l 


Y   X.M  In  X. 
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This  is  true  because  the  coefficients  of  une  equ  tion 
must  be  equal  or  proportional  to  the  coefficients  of  another. 
It  can  now  be  seen  that  successive  iterations  based  on  estimates 
of  M   until  the  determinants  of  the  coefficients  are  equal  to 
zero  to  finalize  the  value  of  the  coefficient  M. 

The  exact  values  of  A  anci  B  can  then  be  calculated  using 
Cramer's  rule. 

The  output  of  the  program  appears  in  two  forms.   One  form 
is  tne  regular  output  of  the  7090  program  which  states  the 
A  and  B  coefficients  ana  tne  exponent  M  for  each  of  the  three 
lines  of  regression.   The  output  of  a  typical  analysis  is 
shown  on  page  13l|f or  examination.   Note  that  with  the  lines  of 
regression  for  points  within  two  and  one  standard  error,  those 
points  which  were  net  used  in  the  determination  ef  the  line 
are  listed  with  the  results.   This  libt  of  points  is  cumulative 
through  the  two  cases.   The  identifying  information  which 
accompanies  each  aata  point  is  placed  in  the  computer  along 
with  the  data  at  the  choice  of  the  user  of  the  program.   The 
second  form  of  output  is  through  the  SC  4020  Charactron.   The 
cooing  for  this  output  is  contained  in  a  subroutine  package 
named  AMPFOS  and  the  form  of  the  output  is  microfilm  copies 
of  the  actual  lines  of  regression.   Although  these  plots 
have  little  value  from  the  practical  standpoint,  due  to  the 
selection  by  the  program's  coding  of  rather  cumbersome  scales, 
the  plots  do  indicate  the  validity  of  the  curve  fit  and  permit 
selection  of  one  of  the  three  lines  as  .being  moot  representative 
of  the  points.   This  output  form  is  shown  on  page  13.5. 
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Results  of  Data  Reduction 

Figures  2  7  thru  50  which  follow  are  plots  of  Nusselt 
Number  versus  Reynolds  Number  for  individual  geometries  and 
angles.   Curves  for  ©  =  90,  75,  60,  45  and  thirty  degrees  are 
plotted  using  experimentally  derived  data,  other  curves  for 
each  geometry  are  also  plotted  using  the  Colburn  and  Grimison 
eguations. 
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APPENDIX  H 

Analysis  of  Frictional  Resistance  Results 

( 
Figures  51  thru  66  are  plots  of  the  computed  friction 

factor  versus  Reynolds  Number  for  individual  geometries  and 

angles  of  inclination.   The  uncertainty  of  data  for  Reynolds 

Numbers  lower  than  2500  was  such  that  no  points  were  used  in 

the  analysis.   The  process  of  curve  fitting   used  was  the 

method  of  least  squared. 

Least  Squares  Analysis 

The  frictional  resistance  data  is  fitted,  to  an  equation 
of  the  form: 

pvpo 

FRIFAC  =  CCNST  (CNRE) 
using  the  method  of  least  squares.   The  logarithmic  equivalent 

of  the  equation  is  taken: 

log  FRIFAC  =  log  CONST  +  EXPO  log  CNRE 
and  the  normal  equations  developed  by  the  usual  method  of 
multiplying  by  the  coefficients  of  the  variables  CONST  and  EXPO 
in  each  case  for  each  data  point  and  summing.   The  normal 
equations  so  developed  are: 

log  FRIFAC  =  (log  CNRE)  EXPO  +  THUMB  log  CONST 

(log  FRIFAC  log  CNRE)  =  (log  CNRE)2  EXPO  +  (log  CNRE)  log  COI^ST 
These  normal  equations  can  be  solved  simultaneously  for  the  log 
of  CONST  and  EXPO: 

EXPO  =  (THUMB)  (log  FRIFAC)  (log  CNRE)-(loq  CNRE)  (log  FRIFAC) 

(THUMB)  (log  CNRE)2  +  (log  CNRE)2 

XCON  =  (log  FRIFAC)  -  ( log  CNRE)  (EXPO) 

(THUMB)  ' 

Then:  C6NST  =  antilog  XCON 

The  parameters  necessary  for  expression  of  friction  factor  as  a 

function  of  Reynolds  Numbers  are  thus  developed. 
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RUNS    3601  -3626 
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Frictional  Resistance  Analysis  Computer  Program ^Description 

The  computer  program  described  herein  is  a  FORTRAN  II 
program  compatible  with  tne  IBM  1620  data  processing  system. 
The  function  of  the  program  is  tc  fit  a  curve  whereby  friction 
factor  is  expressed  as  a  power  of  Reynolds  Number  using  the 
method  of  least  squares.   Ivhile  variable  names,  input,  output 
and  format  statements  are  particularized  to  the  problem  at  hand 
the  methou  is  completely  general  and  can  be  used  in  any 
application  where  a  logarithmic  least  squares  analysis  is 
desired  with  only  minor  program  moaif ications.   It  would  be 
possible  to  include  tnis  program  in  the  basic  data  reduction 
program  to  permit  reduction  and  analysis  in  one  computation. 

Overall  control  of  the  program  is  by  a  DO  loop  indexed  by 
the  number  of  sets  of  data  being  analyzed.   For  each  uata  set 
another  DO  loop  indexed  by  the  number  of  data  points  in  the 
data  set  controls  the  calculations.   For  each  data  set  the 
program  first  sets  the  summations  equal  to  NUME ,  then  accents 
as  input  the  friction  factor  and  Reynolds  nunider  for  the  first 
data  point.   The  natural  logarithms  of  the  input  quantities, 
the  product  of  the  two  logarithms  and  the  square  of  the 
logarithm  of  Reynolds  Number  are  computed.   The  summations 
necessary  to  determine  the  coefficients  of  tne  normal  equations 
are  computed  with  variable  names  being  changed  as  necessary 
for  proper  program  operation.   The  program  tnen  loops  to  repeat 
the  calculations  described  adove  for  each  data  point  until  all 
data  points  in  the  set  hove  be<  n  calculated.   Tne  values  of 
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tne  exponent,  EXPO,  anu  the  natural  logarithm  of  the  coefficient, 

XCON  can  then  be  computed  for  the  particular  data  set.   As  no 

antilogarithm  function  is  available  in  the  FORTRAN  II  library 

the  value  of  the  coefficient,  CONST,  is  calculated  by: 

CONST  =  ,XC0N 
/.hen  tne  values  of  the  coefficient  and  exponent  have  been 

determined  an  output  consisting  of  these  guantities  plus  an 

identifying  coue,  and  the  logarithm  of  the  coefficient  is 

printed. 

The  guantities  which  have  been  computed  so  i ar  represent  the 
parameters  necessary  for  expressing  friction  factor  as  a  function 
of  Reynolds  Number  for  a  particular  geometry  ana  angle  of 
inclination.   In  order  to  correlate  the  values  calculated  with 
the  eguation  developed  using  these  parameters  with  test  results, 
values  of  Reynolus  Number  from  three  thousand  to  eight  thousand 
are  inserted  into  the  eguation  and  the  corresponding  friction 
factor  computed.   The  results  are  printed. 

When  the  calculation  of  the  friction  factors  is  completed 
the  program  loops  to  perform  the  calculations  for  tne  next  data 
set  until  the  limiting  index  is  reached. 

A  program  block  diagram,  listing,  sample  output,  and  output 
summary  follow. 

Computations  using  this  program  were  made  at  C.W.  Post 
College;  data  was  prepared  at  Webb  Institute. 
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C       FORTRAN  TWO  PROGRAM  TO  CORRELATE  FRACTIONAL  RESISTANCE  DATA  BY 
C      LEAST  SQUARES 

1  FORMAT  (U) 

2  FORMAT(  214) 

3  FORMAT ( F9.6,2X,F9.2 ) 

4  FORMAT ( lhlf2X, 57HRUN  CQUE      EXPONENT        CONSTANT      LOGtCONS 
P     IT ANT)     ) 

5  FOkMAT(lh<>,4XtI2,7X,F10.6,5X,F10.6,5X,F10.6) 
9  F0RMAT(1H0,6X,FB.2»13X,F9.6) 

10  FORMAT (  1H0,2X, 37HREYN0LDS  NUMbER       FRICTION  FACTOR  ) 
READ  1,KASES 
DO  7  MA=1,KASES 
SUMFLN=0.0 
SUMREL=O.U 
SMSRG=0.0 
SMFGRE=0.0 
READ  2,NUMB,RNCD 
DO  6  NU=1,NUMB 
READ  3,FRIFAC,CNRE 
FLO=LOGF(FRIFAC) 
RELG=LOGF(CNRE) 
FLGREL=FLG*RELG 
SRELG*RELG**2 
SUMFLG=SUMFLN+FLG 
SUMFLM=SUMFLG 
SMRELG=SUMREL+RELG 
SUMREL=SMRELG 
SMSREL=SMSRG+SRELG 
SMSRG=SMSREL 
SMFLRL=SMFGRE+FLGREL 
SMFGRE=SMFLRL 

6  CONTINUE 
THUMB=NUMB 

EXPO=  (  THUMB*SMFGRE-SUMREl*SUMFLN  )  /  (  THUMb*SMSRG-  (  SUMREL**2  )  ) 
XCON=(  ( SUMFLN-(SUMREL*EXPU)  } /THUMB) 
CPNST=2.71828**XC0N 
PRINT  4 

PRINT  5, RNCDt EXPO, CONST tXCUN 
PRINT  10 
DO  8  NRE=3,8 
DNRE=NRE 

ENRE=DNRE*( 1.0E3) 
FRICT=CONST*(  ENRE**EXPOJ 
PRINT  9,ENPE,FRICT 
8  CONTINUE 

7  CONTINUE 
STOP 
END 


RUN  CODE 

EXPONENT 

CONSTANT 

LOGICUNST/ 

43 

• 

060958 

.033258 

-3.403455 ' 

REYNOLDS  NUMBER 

=RICTION  FACTOR 

3000.00 

.034187 

A  000. 00 

.U55145 

5000.00 

.055901 

61)00.00 

.056526 

7000.00 

.057059 

8000.00 

.057526 

STOP 
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Summary   of    Frictional    Resistance   Analysis 

Run   Code  Exponent  Constant 

14        s  -0.043053  0.065543 

16  0.196590  0.045465 

17  -0.045756  0.445318 
19  -0.079225  0.214984 

24  -0.024007  0.100758 

26  -0.119514  0.037337 

29  -0.081586  0.261658 

33  0.057896  0.040784 

34  -0.089713  0.270171 
36  -0.028274  0.137797 
39  -0.086865  0.U31623 

43  0.060968  .          0.033258 

44  0.210122  C. 017001 

45  0.053331  0.074656 
49  -0.055798  0.297399 
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Run  Code             14           16           17  19 

Reynolds  Friction  Friction     Friction  Friction 

Number  Factor  Factor       Factor  Factor 

3000  0.046433  0.219407  0.308722  0.114006 

4000  0.045861  0.^32174  0.304685  0J11437 

5000  0.045422  0.242585  0.301590  0.109484 

6000  0.045067  0.251438  0.299084  0.107914 

7000  C. 044769  0.259174  0.296982  0.106604 

8000  0.044513  0.266068  0.295173  0.105482 

Run  Code  24  26  29 

Reynolds  Friction  Friction      Friction 

Number  Factor  Factor        Factor 

3000  0.083138  0.097209  0.136263 

4000  0.082566  0.110610  0.133102 

5000  0.082125  0.103329  0.130701 

6000  0.081766  0.105605  0.128771 

7000  0.081464  0.107569  0.127162 

8000  0.081204  0.109299  0.125784 

Run  Code             33           34           36  39 

Reynolds  Friction  Friction     Friction  Friction 

Number  Factor  Factor        Factor  Factor 

3000  0.064834  0.131732  0.149752  0.215310 

4000  0.065923  0.128376  0.148539  0.209996 

5000  0.066780  0.125831  0.147605  0.205965 

6000  0.067489  0.123790  0.146846  0.202728 

7000  0.068094  0.122090  0.146207  0.200032 

8000  0.068622  0.120636  0.145656  0.197725 
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Run  Code  43  44  46  49 

Reynolus  Friction  Friction  Friction  Friction 

Number  Factor  Factor  Factor  Factor 

3000  0.054187  0.091434  0.114421  0.190248 

4000  0.055145  0.097132  0.116190  0.187219 

5000  0.055901  0.101795  0.117581  0.184902 

6000  0.056526  0.105770  0.118730  0.183031 

7000  0.051059  0.109252  0.119710  0.181463 

8000  0.057526  0.112361  0.120568  0. 180116 
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